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(54) Method for determining a partition in order to insert a watermark, and associated insertion 
and decoding methods 



(57) A method of determining a partition into distinct 
regions of a set of coefficients representative of at least 
a part of an image in order to insert a watermark com- 
prises the following steps: 

generation (E501 ) of a centred pseudo-random se- 
quence (w) equal in size to the cardinal number (N) 
of the said set of coefficients, formed from centred 
pseudo-random sub-sequences; 
modulation (E502) of the said set of coefficients by 
the said centred pseudo-random sequence (w) in 
order to insert the same information bit on the said 
set of coefficients; and 

checking (E504-E508), over each region : of a crite- 
rion of detectability of the information bit inserted by 
modulation. 

Use especially for determining an adaptalive parti- 
tion as a function of a criterion of detectability of an in- 
formation bit inserted on each region. 
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Description 

[0001] The present invention concerns a method of determining a partitioning of a digital image for inserting a wa- 
termarking signal, a method of inserting a watermark, a method of processing an image and an associated method of 
decoding a watermarking signal. 

[0002] Correlatively, it concerns a device for determining a partitioning of a digital image for inserting a watermarking 
signal, a device for inserting a watermark, a device for processing an image and a device for decoding a watermarking 
signal. 

[0003] The present invention lies in general terms in the technical field of watermarking of digital images, and more 
particularly still images. 

[0004] Watermarking digital data makes it possible to protect these data, for example by associating copyright infor- 
mation with them. 

[0005] In its general principe, watermarking consists of inserting an indelible watermark in digital data, which can be 
likened to encoding additional information in the data. 

[0006] Decoding this additional information makes it possible to check the additional information which has been 
inserted. 

[0007] This inserted watermark must consequently be both imperceptible, robust to certain distortions applied to the 
digital image and reliable to detect. 

[0008] Conventionally, a usual technique of inserting a watermarking signal in a digital image consists of using a 
linear modulation model in which at least one subset of coefficients representing the digital image is modulated ac- 
cording to this linear model using a weighting coefficient. 

[0009] By denoting as X = {Xj, 1 < i < N) a set of the coefficients representative at least of part of a digital image, and 
as w = {w j; 1 < j < P} a watermark of size P < N, a pseudo-random signal with known distribution and average of zero, 
the linear insertion formula is : 



XV = ^ + 5^ Wj, with 1 < j < P, 

in which {X y 1 < j < P} is a subset of the set of coefficients X, b is an information bit, and ou is a weighting coefficient, 
also called modulation amplitude. 

[0010] The detection of the watermark then consists in detecting whether or not the pseudo-random sequence w 
has been inserted into a set of coefficients. This detection is carried out without using the original image and can be 
based on a standardised statistical test which makes it possible to calculate a probability of correct detection. 
[0011] Such an insertion technique makes it possible, by insertion of a watermark, to insert a single information bit 
since the response from the detector is binary (yes/no). 

[0012] In order to insert a larger number of information bits into the digital image, in particular when a code of Q bits 
is wanted, indicating, for example, the name or the address of the owner or of the author of the image, it is necessary 
to reiterate the previously described insertion method as many times as there are information bits to be inserted. 
Typically, in order to insert a binary signal, either b = 1 , or b = -1 is used. 

[001 3] Put another way, Q subsets of coefficients have to be chosen and the modulation of these subsets has to be 
carried out by choosing Q watermarks. 

[0014] Separate subsets of coefficients are preferably chosen such that the modulation operations are not superim- 
posed on one another, which could disturb the detection or cause troublesome visual effects. 

[0015] It is a matter, consequently, of choosing a partition of the coefficients representative of the digital image into 
Q separate subsets, each carrying one information bit. 

[0016] Numerous known methods use a technique of inserting a watermarking signal of given size by spectrum 
spreading. The disadvantage of these methods, describe for example in the article entitled "Secure spread spectrum 
watermarking for multimedia" by I.J. COX el al, in Proc. ICIP, pages 243-246, Sept. 1996 and in the article entitled 
"Digital watermarking of raw and compressed video" by F. HARTUNG et al, in Proc. SPIE 2952 : Digital Compression 
Technologies and Systems for Video Communication, pages 205-213, October 1996, is that they use an arbitrary 
partitioning of the image into blocks of fixed size without any guarantee of detectability of the modulations effected on 
each of the blocks. 

[0017] The aim of the present invention is to propose a method of determining a partitioning of a digital image which 
is adapted to the image, ensuring detectability of the bits inserted for a technique of insertion by modulation. 
[0018] The present invention thus relates to a method of determining a partitioning into distinct regions of a digital 
image in order to insert a watermarking signal by modulation of coefficients of said regbns representing a digital image. 
[0019] In accordance with the invention, this partitioning into regions is effected by an adaptive partitioning as a 
function of a criterion of detectability of an information bit inserted in each region. 
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[0020] In this way an adaptive partitioning is performed which has the advantage of adapting to the content of the 
image signal and guaranteeing detectability of the inserted watermarking signal, unlike the choice of an arbitrary par- 
titioning into fixed blocks. 

[0021] This spatial adaptivity is necessary so as the content of the image is not spatially homogenous. 
5 [0022] This is because, for a given modulation amplitude, the greater the length of the sequence of the watermarking 
signal inserted, and consequently the greater the size of the region and therefore the number of modulatable coeffi- 
cients, the more reliable is the detection. 

[0023] It can in fact be shown that, when a correlation calculation is used for effecting the detection, the response 
of the detector is proportional to P 1/2 where P is the length of the sequence of the watermarking signal. 
io [0024] Moreover, the response of the detector is also dependent on the image signal itself: The greater the variance 
in the coefficients representing the image, the weaker the response of the detector. 

[0025] In conclusion, for a given modulation amplitude and image, there exists a minimum length of the watermarking 
signal necessary for ensuring a given detection probability, this minimum length also depending on the local content 
of the image signal in the region in which the watermarking signal is inserted. 
75 [0026] The adaptive partitioning effected by the method according to the invention guarantees the detectability of 
the watermarking signal inserted in the different regions determined by this partitioning. 

[0027] Preferably, the detectability criterion is calculated using a maximum value of the weighting coefficient ensuring 
imperceptibility of the watermarking signal for the modulation of the set of coefficients representing said digital image. 
[0028] This is because, the greater the modulation amplitude, that is to say the value of the weighting coefficient, 
20 the more reliable is the detection in that a better probability of delection can be obtained. Nevertheless, the greater 
this modulation amplitude, the more this modulation may be visible in the watermarked image. 

[0029] It is therefore judicious to choose a modulation amplitude equal to a maximum amplitude, generally denoted 

JND (Just'Noticeable Difference), beyond which the modulation becomes visibly perceptible in the image. 

[0030]". It can be shown that this maximum amplitude also depends on the length of the sequence of the watermarking 

25 signal and that it decreases with this length. 

[0031] By choosing a maximum value of the weighting coefficient ensuring imperceptibility of the watermarking signal 
for the modulation of the set of coefficients representing said digital image, the value of the weighting coefficient on 
each region of the partitioning, and with a size less than the cardinal number of all the coefficients representing the 
image, will be much less than the maximum modulation amplitude which guarantees imperceptibility of the modulation 

30 on this region. 

[0032] The partitioning thus performed also ensures invisibility of the modulation effected on each of the regions 
obtained by the partitioning. 

[0033] According to an advantageous characteristic of the invention, the detectability criterion is calculated using, 
for each coefficient to be modulated, a weighting coefficient according to a law of the form cxj = k^, where kj is a 
35 modulation factor dependent on the coefficients situated in a vicinity of the coefficient considered on said region and 
a v is equal to said maximum value of the weighting coefficient. 

[0034] It is thus possible to exploit the possibilities of visual "masking" of the modulation by the image signal itself. 
Typically, use is made of the fact that, the greater the activity of the image signal , in the sense for example of a significant 
local variance in the coefficients, the less visible will be the modulation. 
40 [0035] It is thus possible to locally increase the weighting coefficient and therefore the modulation amplitude, to the 
benefit of detectability, without impairing the imperceptibility of the inserted modulation. 

[0036] According to another preferred characteristic of the invention, the method of determining a partitioning includes 
a step of applying a distortion before calculating the detectability criterion. 

[0037] According to another preferred characteristic of the invention, the method of determining a partitioning includes 
45 a prior step of comparing the size of each region of the partitioning with a minimum size corresponding to the minimum 

size of a statistically significant sample for detecting a watermarking signal inserted in said region. 

[0038] indeed, because of the calcu lations of a statistical nature performed when the watermarking signal is detected, 

notably when standardised hypothesis tests are used, it is necessary to consider a minimum boundary for the size of 

each region ol the partitioning guaranteeing the validity of the detection tests. 
so [0039] According to one particularly practical embodiment of the invention, for each region of the partitioning, all the 

coefficients of said region are modulated by inserting the watermarking signal, a detectability amplitude is calculated 

from said detectability criterion and said detectability criterion is validated by comparing said detectability amplitude 

with a predetermined threshold value. 

[0040] The detectability amplitude can be the result of a calculation of correlation between the set of modulated 
55 coefficients and the watermarking signal or the result of a standardised statistical test such as the one described in 
the article entitled "A method for signature casting on digital images" by I. PITAS, in Proc. ICIP, pages 215-218, Sep- 
tember 1996. 

[0041] Practically, the detectability criterion can be validated if a minimum modulation length calculated for ensuring 
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detectabilily of the inserted watermarking signal in less than or equal to the cardinal number of said region 
[0042] According to a preferred embodiment of the invention, the method comprises the following steps: 

- generation of a centred pseudo-random sequence equal in size to the cardinal number of a set of coefficients 
representative of at least a part of the image, formed from centred pseudo-random sub-sequences 

- modulation of the said set of coefficients by the said centred pseudo-random sequence in order to insert the same 
information bit on the said set of coefficients; and 

- checking, over each region, of a criterion of detectabilily of the information bit inserted by modulation. 

[0043] It is thus possible to achieve a partition of a set of coefficients into distinct regions, by carrying out an adaptive 
partitioning on the basis of a criterion of detectabilily of an information bit inserted on each region 
[0044] Such an adaptive partitioning exhibits the advantage of adapting to the content of the image signal and of 
guaranteeing the detectabilily of the watermark inserted on the different regions determined by this partition in contrast 
to the choice of a random partitioning. 

[004S] In practice, for each region of the partition, a crilerion of detectabilily of an information bit inserted by modu- 
lation of the coeft.cients of the said region is checked. This criterion may consist in the comparison of a result of a 
standardised statistical test with a threshold value corresponding to a probability of detection fixed in advance 
[0046] By virtue of the method in accordance with the invention, it is possible to generate a centred pseudo-random 
sequence on the basis of the merging of centred pseudo-random sequences of smaller sizes, and to modulate all the 
available coefficients in order to achieve the insertion of a watermark. 

[0047] By modulating all the coefficients with the same information bit and a centred pseudo-random sequence which 
is broken down into pseudo-random sub-sequences which are also centred on each region of the partition produced 
over the set of coefficients, the steps of generating a pseudo-random sequence and of modulation do not have to be 
reiterated at each new region tested. 

[0048] According to one preferred characteristic of the invention, the method of determining a partition comprises a 
criterion^ ^ aPP ' ym9 ® distortion to the said set of m °dulated coefficients, before the step of checking of a detectability 

[0049] It is known to apply a predetermined distortion to a set of modulated coefficients and to carry out the check 
on the detectabilily criterion after this pre<fistortion, in such a way that it is possible to guarantee the detection of the 
watermark inserted on a partition of the coefficients for a predetermined distortion. 

[0050] By virtue of the invention, a pre-distortion can be applied to the set of modulated coefficients, in a unique way, 
and does not need to be reiterated at each new region tested in the partition. 

[0051] In one particularly practical embodiment of the invention, the set of coefficients is a set of spectral coefficients 
originating from a spectral transformation of a digital image, and the step of applying a distortion includes the following 
sub-steps: M 

- inverse spectral transformation of the said set of modulated coefficients in order to obtain a watermarked imaqe- 
application of a distortion to the said watermarked image; and ' 

- spectral transformation of the watermarked image in order to obtain a set of modulated coefficients after distortion. 

[0052] The saving obtained in terms of complexity of calculation is all greater in this case whenever the application 
of a distortion requires several steps of processing of the set of coefficients, and particularly the application of an 
inverse spectral transformation in order to obtain the watermarked image and to apply a distortion, such as compres- 
sion/decompression to it, and then the application of a forward spectral transformation of the image again so as to 
check the criterion of detectability in the spectral domain. 

[0053] According to another particularly practical characteristic of the invention, the centred pseudo-random se- 
quence is formed from pseudo-random sub-sequences of identical size. 

[0054] Such a pseudo-random sequence is particularly well adapted to the application of a bottom-up or top-down 
iterative partitioning of the quadtree type which makes it possible iteratively to combine regions of identical sizes or 
iteratively to divide regions into sub-regions of identical sizes. 

[0055] For preference, the pseudo-random sub-sequences have a size greater than or equal to a minimum number 
of coeff.c.ents forming a region of the said partition of a size which is statistically significant for validating the said 
criterion of detectability over the said region. 

[0056] It is thus possible theoretically to form a maximum partition of the set of coefficients with the same unique 
centred pseudo-random sequence, each pseudo-random sub-sequence being centred on the regions with the smallest 
size which is statistically admissible for inserting a watermark by modulation into the set of coefficients 
[0057] According to one practical embodiment, which is inexpensive in terms of calculations, at the generation step 
each centred pseudo-random sub-sequence is created by generating half of the pseudo-random samples of the said 
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sub-sequence via a pseudo-random number generator and by symmetrising the said pseudo-random samples gener- 
ated. 

[0058] For preference, in this embodiment, at the generation step, each centred pseudo-random sub-sequence is 
created by moreover permutating the said pseudo-random samples using a secret key. 
s [0059] This permutation operation makes it possible better to preserve the intrinsic random-sequence nature of a 
sub-sequence created by symmetrising pseudo-random samples, and not to diminish the imperceptibility of the inserted 
watermark. 

[0060] According to another practical embodiment of the invention, at the generation step, each centred pseudo- 
random sub-sequence is created by generating pseudo-random random samples via a pseudo-random number gen- 

io erator and by redistributing over each sample the sum of the said pseudo-random samples generated. 

[0061] According to a preferred characteristic of the invention, a partitioning by iterative division of the digital image 
is performed and : for each region at a rank of said partitioning, a partitioning of a directly higher rank of said region is 
validated if and only if there exist at least two subregions of said region for which said detectability criterion is validated. 
[0062] Thus, when it is attempted to divide a region into subregions, this partitioning of a higher rank is validated 

J5 only if the number of regions in which a watermarking signal can actually be inserted whilst complying with the detecta- 
bility criterion is actually increased. 

[0063] Otherwise, it is preferable to keep the region before division, of greater size, in order to increase the detecta- 
bility of the watermarking signal. 

[0064] According to a preferred embodiment of the invention, for each region at a rank of said partitioning, the value 
20 of a capacity of the digital image is incremented when the partitioning of directly higher rank is not validated, the said 
capacity being equal to the cardinal number of the distinct regions obtained by the said partitioning on which the de- 
tectability criterion is validated. 

[0065] It is thus possible to calculate the capacity of the digital image which thus corresponds to the number of 
information bits which can be inserted in the digital image. 

25 [0066] By virtue of the adaptive partitioning performed by the invention, by reiterating the partitioning of each region 
as long as the detectability criterion is complied with, there is obtained the maximum capacity of the digital image which 
corresponds to the watermarking capacity of this image by the technique of insertion by modulation. This watermarking 
capacity can be defined as the maximum number of information bits which can be inserted imperceptibly in the image 
with a guarantee of detectability on decoding. 

30 [0067] According to another embodiment of the invention, for each region at a rank of said partitioning, the value of 
a capacity is incremented when the partitioning of directly higher rank is not validated, said value of the capacity is 
compared with a fixed capacity value and the partitioning is reiterated for a region having a detectability amplitude 
which is the highest amongst all the other regions to be processed when said capacity value is less than said fixed 
capacity value. 

35 [0068] This partitioning method makes it possible, in this embodiment, to find the best possible partitioning in terms 
of detectability when a predetermined number of information bits have to be inserted, equal to said fixed capacity value. 
[0069] By selecting, during the partitioning, the regions for which the detection amplitude is the highest, the present 
invention makes it possible to determine the partitioning which, a priori, guarantees the best detectability on decoding 
of the watermarking signals making it possible to insert a predetermined number of information bits. 

40 [0070] According to another embodiment of the invention, a partitioning is effected by iterative merging of the digital 
image and, for each region of said partitioning for which said detectability criterion is not validated, said region is merged 
with another region of said partitioning. 

[0071] Preferably, in this embodiment, said other region of the partitioning is if possible a region for which said de- 
tectability criterion is not validated. 
45 [0072] It is thus possible to increase the capacity of the digital image by joining the regions in which the detectability 
criterion is not satisfied. 

[0073] According to another aspect of the invention, a method of inserting, into a digital image, a watermark com- 
prising one or more information bits, comprises a method of determining a partition in accordance with the first aspect 
of the invention, and a step of inserting the information bit or bits by modulation of the coefficients of one or more 
50 regions of the said partition. 

[0074] This insertion method exhibits advantages similar to those of the method of determining a partition which it 
implements, and in particular makes it possible to insert a watermark in an adaptive and reliable way as regards its 
detection. 

[0075] The present invention also concerns a method of decoding a watermarking signal in a digital image, inserted 
55 by modulation of coefficients representing said image in distinct regions forming a partitioning of said digital image. 
[0076] In accordance with the invention, this decoding method includes a step of determining a partitioning into 
distinct regions of the digital image to be decoded by effecting an adaptive partitioning as a function of a criterion of 
detectability of an information bit obtained by demodulation of the coefficients on each region. 
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[0077] Thus this decoding method makes it possible to find, at the decoder, the adaptive partitioning performed for 
inserting a watermarking signal. 

[0078] This is because, when the partitioning has been performed in an adaptive fashion as described above by a 
method of determining a partitioning in accordance with the first aspect of the invention, this partitioning, which depends 
on the content of the image, is not known to the decoder. 

[0079] The decoding method according to the invention then consists of testing all the acceptable partitionings and 
verifying the detection of a watermarking signal on each of the regions associated with this partitioning. 
[0080] According to a preferred characteristic of the invention, the decoding method includes a prior step of comparing 
the size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistica ly 
significant sample for the detection of a watermarking signal inserted in said region. 

[0081] According to a practical embodiment of the decoding method according to the invention, for each region of 
the partitioning, a detection amplitude is calculated from said detectability criterion and said detectability criterion is 
validated by comparing said detection amplitude with a threshold value predetermined for decoding 
[0082] The decoding method thus makes it possible to find the adaptive partitioning effected when the watermarking 
signal is inserted. 

[0083] According to advantageous characteristic of this embodiment in the invention, the threshold value predeter- 
mined for decoding is less than said predetermined threshold value used during the method of determining a partition- 
ing. 

[0084] This is because, at decoding, the detection test is performed on coefficients which have a priori been modu- 
lated and possibly made noisy. These coefficients are different from the original coefficients, and the resull of the 
detectability test at decoding will be different from that performed in order to determine the partitioning. By choosing 
different thresholds on decoding and coding, it is possible to take into account the tendency for variation in detection 
amplitude due to an increased variance in the coefficients which have been modulated and possibly made noisy 
[0085] Correlatively, the present invention also concerns a device for determining a partitioning into distinct regions 
of a digital image for inserting a watermarking signal by modulating coefficients of said regions, comprising partitioning 
means adapted to effect an adaptative partitioning as a function of a criterion of detectability of an information bit 
inserted on each region. 

[0086] According to a preferred embodiment, this device for determining a partition includes: 

means for generating a centred pseudo-random sequence equal in size to the cardinal number of the said set of 
coefficients, formed by centred pseudo-random sub-sequences; 

means for modulation of a set of coefficients by the said centred pseudo-random sequence in order to insert the 
same information bit on the said set of coefficients: and 

means for checking, over each region, a criterion of detectability of the information bit inserted by modulation. 

[0087] The present invention also envisages a device for inserting into a digital image a watermark comprising one 
or more information bits, comprising a device for determining a partition in accordance with the invention, and means 
for inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 
[0088] Correlatively, the present invention also concerns a decoding device adapted to implement the decoding 
method according to the invention. 

[0089] The present invention also concerns a method of processing an image, comprising the following steps: 
inputting an image; 

inserting a first watermarking information into said inputted image; 
applying a distortion to the image obtained in said insertion step; 

detecting said first watermarking information from the image obtained in said distortion applying step; and 
determining an amount of a second watermarking information which can be inserted into said inputted image in 
accordance with the result of the detection obtained in said detection step. 

[0090] It further comprises the step of displaying the amount of the watermarking information determined in said 
determining step. 

[0091] Correlatively, the present invention also concerns an image processing apparatus, comprising: 
means for inputting an image; 

means for inserting a first watermarking information into said inputted image; 
means for applying a distortion to the image obtained by said insertion means; 

means for detecting said first watermarking information from the image obtained by said distortion applying means; 
and 
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means for determining an amount of a second watermarking information which can be inserted into said inputted 
image in accordance with the result of the detection obtained by said detection means. 

[0092] This apparatus further comprises the means for displaying the amount of the watermarking information de- 
5 termined by said determining means. 

[0093] The present invention also concerns a method of inserting a watermarking signal which can take into account 
at the same time predetermined imperceptibility and detectability constraints in order to ensure invisibility of the wa- 
termarking signal in the image and correct detection at the time of decoding. 

[0094] In the article entitled "A method for signature casting on digital images" by I. PITAS, in Proc. ICIR pages 
io 215-218, September 1996, there is described a method of inserting a watermarking signal in a digital image in which 
the weighting coefficient a is determined so as to ensure detectability of the signal with a fixed probability 
[0095] This weighting coefficient must be greater than a minimum value, which can be referred to as the detection 
amplitude, so as to allow the detection of the inserted signal with the correct detection probability level. 
[0096] This minimum value or detection amplitude depends notably on the length of the watermarking signal and on 
15 the required detection probability level. 

[0097] However, the weighting coefficient is fixed in the aforementioned document without taking into account the 
concept of psychovisual amplitude which guarantees invisibility of the watermarking signal inserted in the digital image. 
[0098] To this end the present invention relates to a method of inserting a watermarking signal in a set of coefficients 
representing a digital image, in which at least one subset of coefficients is modulated according to a linear model using 
20 a weighting coefficient, the weighting coefficient being greater than a minimum value determined as a function notably 
of the length of the watermarking signal and a detection probability level. 

[0099] According to the invention, this insertion method includes a step of calculating a maximum value of the weight- 
ing coefficient according to the length of the watermarking signal ensuring imperceptibility of the watermarking signal. 
[0100] Calculating a maximum value for the weighting coefficient, which can be referred to as visual amplitude, 
25 makes it possible to take into account, when a watermarking signal is inserted in a set of coefficients representing a 
digital image, constraints related to the imperceptibility of the inserted watermarking signal. 

[0101] Such an insertion method makes it possible to determine the conditions which make possible the insertion of 
a watermarking signal in an image whilst meeting predetermined imperceptibility and detectability criterion or on the 
other hand to decide if, under certain conditions, the insertion of the watermarking signal meets the predetermined 

30 imperceptibility and detectability criterion. 

[0102] According to a practical characteristic of implementation of the present invention, the insertion method also 
includes a step of determining compatible values of the length of the watermarking signal, of the weighting coefficient 
and of the probability level so that the value of the weighting coefficient determined is less than or equal to said maximum 
value and greater than or equal to said minimum value for said determined values of the length of the watermarking 

35 signal and of the detection probability level. 

[0103] By satisfying the aforementioned inequality for the weighting coefficient, the conditions of detectability and 
imperceptibility are satisfied at the same time. 

[0104] It is thus possible to determine, in a correlated fashion, the value of the weighting coefficient, of the length of 
the watermarking signal and of the detection probability level in order to satisfy this inequality 
40 [01 05] According to a preferred characteristic of the invention, during the step of calculating the said maximum value, 
a function independent of said set of coefficients to be watermarked is used. 

[0106] Thus the model used for calculating the maximum value or the visual amplitude is independent of the image 
to be watermarked. 

[0107] In a preferred embodiment of the invention, the coefficients are spatio-frequency coefficients obtained by 
45 spatio-frequency transformation of this image, the calculation of said maximum value of the weighting coefficient de- 
pending on the type of transformation used. 

[0108] It is normal, before the insertion proper of a watermarking signal to apply a spatio-frequency transformation 
to the image, also referred to as transformation into sub-bands of the image, of the discrete cosine transformation 
(DCT) or wavelet transformation type, so that the modulation for inserting the watermarking signal is implemented on 
50 spatio-frequency coefficients of a sub-band of the image. 

[0109] The insertion method according to the invention applies particularly well to the conventional techniques of 
insertion by modulation inserted in a transformed domain obtained by one ol the aforementioned spatio-frequency 
transformations. 

[0110] In practice, said calculation function used for calculating the maximum value or visual amplitude c^asa 
55 function of the length P of the watermarking signal w is of the type: 
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p p {£[iivr]) [ 



where 



a base (S) is a base value, dependent on the transformation used S and on the (base) sub-band considered for the 
insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of said sub-band 
is modulated, 
10 (3 is strictly greater than 2, and 

E[lwlP] is the mathematical expectation of the function IwlP 



[0111] According to another preferred characteristic of the invention, said minimum value is determined by using a 
calculation function also depending on the variance of the coefficients to be modulated, on the distribution of the wa- 
15 termarking signal and possibly on the variance of an additive noise decorrelated from the coefficients modelling any 
distortions applied to coefficients. 

[0112] Thus the minimum value of the weighting coefficient or detection amplitude can be calculated by taking into 
account notably the distortions which the image may undergo, notably during normal processings consisting of com- 
pression of digital images, in so far as such distortions can be modelled in the form of an additive noise decorrelated 
20 from the coefficients representing the image. 

[0113] In practice, the said function of calculating the detection amplitude a D as a function of the length P of the 
watermarking signal w is of the type 



25 



where a£ is equal to the variance of the coefficients to be modulated, is equal to the variance of the additive 
noise and a, b and c are constants which depend on the distribution of the watermarking signal and on the required 
30 detection probability level. 

[01 1 4] According to a preferred characteristic of the invention, the watermarking signal is a pseudo-random sequence 
with a predetermined distribution with a zero mean. 

[0115] The use of such a watermarking signal increases the robustness of the inserted watermarking signal with 
respect to an intentional attack aimed at extracting or modifying the watermarking signal, typically by statistical analysis, 
35 the kernel of the pseudo-random sequence produced from a predetermined distribution also being able to be given by 
a secret key. 

[0116] In a particularly advantageous embodiment of the invention, the method of inserting a watermarking signal 
includes the following steps: 

40 - calculating a maximum value of the weighting coefficients for a length of the watermarking signal equal to the 
cardinal number of the set of modulatable coefficients ensuring imperceptibility of the watermarking signal at the 
limit of perceptibility; 

calculating the length of the watermarking signal for a minimal value of the weighting coefficient equal to said 
maximum value calculated and for a predetermined detection probability level; and 
4 5 - comparing said calculated length with the cardinal number of the set of modulatable coefficients. 

[0117] Thus, by determining the maximum value of the weighting coefficient at the limit of perceptibility, that is to say 
in the eventuality of all the coefficients available being modulated, it is possible to determine the minimum length of 
the watermarking signal which makes it possible to obtain a correct detection probability with a weighting coefficient 
50 equal to the maximum value calculated. 

[0118] The comparison of this length thus calculated with the cardinal number of all coefficients available in order to 
effect the insertion makes it possible to determine whether this insertion is actually possible with the required detection 
probability level. 

[0119] In addition, if the calculated length is less than the cardinal number of the set of coefficients, it can be certain 
55 that the imperceptibility criterion will indeed be met by modulating a subset of coefficients with a cardinal number equal 
to this calculated length, since the value of the weighting coefficient was calculated in the eventuality of all the available 
coefficients actually being modulated. 

[0120] In practical terms, if. at said comparison step, said calculated length is less than or equal to the cardinal 
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number of the set ot modulatable coefficients, a watermarking signal with a length at least equal to said. Hdeulated 
.ength is inserted by modulating a subset of coefficients with a cardinal number equal to sa.d length, according to a 
linear model using a weighting coefficient equal to said calculated maximum value^ 

[01211 Conversely, if at said comparison step said calculated length is greater than the card.nal number of the set 
5 of modulatable coefficients, there is calculated the detection probability level obtained for a m.n.mum value of the 
weighting coefficient equal to said maximum value calculated lor a watermarking s.gnal length equal to the card.nal 
niimhpr of all the modulatable coefficients. 

JEST Thus the decision to insert a watermarking signal or not on the set of modulatable coeffic.ents can be taken 
in the theoretical knowledge of the detection probability level obtained by default. 
,o 0 23 in another practical embodiment of the invention, the insertion method includes a step of «^^>£>£ 
Lid va ue of the length of the watermarking signal determined so that said minimum value of the we.ght.ng contents 
Ts equal to said maximum value of the weighting coefficients for a length of the watermarking signal equal to sa.d 
threshold value and a predetermined detection probability level. 

r0124l It is thus possible to determine the minimum length threshold which the watermarking s.gnal must meet ,n 
,5 order to simultaneously meet the two criteria of perceptibility and detection when the watermark.ng s.gnal ,s .nserted 

mi^f'coTreSvely, the present invention also concerns a device for inserting a watermarking signal in a set of 

Senu "presenting a Sigita. image, in which at least one subset o, coefficients is modulated 

model using a weighting coefficient, the weighting coefficient being greater than a m.n.mum value determined as a 

20 function notably of the length of the watermarking signal and ot a detection probab.l.ty level 

SST This insertion device comprises means of calculating a maximum value of the weighting coeff.cent as a f unc- 
ion of the length of the watermarking signal ensuring imperceptibility of the watermarking signal. 
[0127] This insertion device has characteristics and advantages similar to those described prev.ously since it » 
adapted to implement the insertion method according to the invention. 

2 s [0128] The present invention also concerns a computer, a digital image process.ng apparatus, a digital p inter, a 
Sigita photographic apparatus and a digital camera adapted to implement the method of determining a 
the methods of inserting a watermark, the method of processing an image and the decod.ng method ,n accordance 

"These^evices for determining a partitioning, inserting a watermark, processing an image and deeding this 
computer thisdigita. image processing apparatus, this digital printer, this digital photographic apparatus and this d.g.tal 
SmeS have characteriLs and advantages similar to those described with reference to the methods wh.ch they 

[0?^ The present invention also relates to a computer program product and to an ^ m ^ Xox ^^^^ 
can be read by a computer or a microprocessor, possibly removable, which stores a program .mplement.ng the method 
o?determ?ning a partitioning and/or inserting a watermark and/or processing an image and/or decod.ng ,n accordance 

ToiVlT ^^particularities and advantages of the invention will also emerge from the following description of the 
invention. 

[01 32] In the accompanying drawings, given by way of non-limitative examples: 

- Figure 1 is an overall diagram illustrating a device for inserting a watermark into a digital image >> 
. Figure 2 schematically illustrates a step o, spectral decomposition employed when mert.ng a 

- Figure 3 is an algorithm of the method of determining a partitioning in accordance wrth a f.rst embod.ment of the 
invention; 

45 - Figure 4 is an example ot partitioning of a region by division; h om hnHim»m nt 

- Figure 5 is an algorithm of the method of determining a partitioning in accordance wrth a second embod.ment of 

- Figure Vis an algorithm of the decoding method according to a first embodiment of the invention; 

- Figure 7 is an algorithm illustrating the application of the method of determining a partrt.on.ng of F.gure 1 to a d.g.tal 

- RgTe 8 is an algorithm illustrating the application of the decoding method of Figure 4 to a digital image, 

- Figure 9 is a block diagram illustrating a variant embodiment of the method of determ.n.ng a partrt.on.ng .n accord- 
ance with a third embodiment of the invention; . 

- Figure 10 illustrates the application of the method of determining a partition to a digital .mage accord.ng to a 

" . K£™ 1 S r5S^?^S£ of determining a partition in accordance with a first way of the preferred 

- igZSZZZZZ oHhe method of determining a partition in accordance with a second way o, the preferred 
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embodiment of the invention ; 

Figure 13 is an algorithm of the insertion method according to one embodiment of the invention- 
Figure 14 illustrates, by means of a graph, the insertion method of the figure 13 and 

Figure 15 ,s a block diagram illustrating a device adapted to implement the methods in accordance with the inven- 



™ J? 1 f deSCription Wi " first of a!l be 9 |ven ' b V reference to Figure 1 , of the insertion of a watermark into a set of 
coefficients representative of a digital image I. 

Hnml J* 1 ' 5 di f 81 T 96 ! °T be re P resen,ed b y a set of coefficients, either in the spatial domain or in a transformed 
domain, the coefficients ,n this latter case being hybrid coefficients, that is to say coefficients located both in the soatfal 
domain and ,n the frequency domain. Such a representation of the image is obtained, for examp.e, by us ng a sS 
r013 5 r 9 H^? h ° m a rt SPaC r r qU6nCy deCOm P° si,ion of the ^e. for example a discrete wavelet decomp'ositon 
ElSi^ T , watermark ^ achieved via a spread-spectrum insertion technique, by modulation of 

dSaNmage ! mentation of the image, obtained by the space-frequency transformation oHhe 

2 01 tS JIT °' ° Xam T • 3 d ' SCrete WaV6,et decom P° sition S can be used, represented diagrammatica.ly ,n Figure 
feite^tedTIwTT f ™T ™ * ^ <° ,h6 im age-processin g field, and on.y the principle thereof s 

reiterated below. It makes it possible to segment the image into sub-bands of frequencies and to obtain hybrid coeffi- 
cient^ that ,s to say spectra, coefficients which are also located in space, here in the plane of the image 

Jk k 'r 9 ? C ° nSiS,S °' 3 Seri6S ° f di9i,a ' Samp ' eS The ima 9 e 1 is represented, for example by a series of 
™ 3 81 Th * 7 7; eSen,,n9 ° ne pixel of the ima 9e I , which may be a black-and-white image wJh 256 grey levels 
[0138] The multi-resolution spectral decomposition means consist of a sub-band decomposLn circuit or analyJs 
c,rcu,t, formed by a set of analysis filters, associated respectively with dividers-by-two. This decomposition circuit fites 
he image s.gna. I ,n two directions, into sub bands of high spatial frequencies and of low freq'uenciesThe 

of reso,ut!r r """^ * ^ ** dOWn the im " e ' int ° • ub * an * accordin 9 to '-els 

!m 4o! , 6Xam M l6, , ,he ima ° e ' iS br ° ken dOWn hSre in, ° sub -bands with a level of decomposition d equal to 3 

Si??f ■» 1 yS,£ l reCe ' VeS ,h6 ima9e Si9nal ' and fM,erS il throu 9 h ,wo respectively low-pass and high-pass 

dig.ta filters, in a first direction, for example horizontal. After passing through dividers-by-two, the resultant filled 

vertiS 7 V T ,i,,ered ^ tW ° reSP6C,iVely '° W - paSS and hi9h - paSS fi,ters ' in a second di -.ion. for LamSe 
vertical. Each signal is aga.n passed through a divider-by-two. Hence four sub-bands LL„ LH„ HL, and HH, are 

22« ?l ° UtPUt ° f thiS ,irSt analVSiS bl0Ck ' Wi,h the hi 9 hest reso,ution in ^e decomposition 1 
ll,h lU^ S ; band ,l; Ll 1 inC i UdeS the low - ,rec ' uenc y components in the two directions of the image signal I The 
second Hii ; ' n r f h S IOW - frec < uenc y components along a first direction and high-frequency components along a 

dSltnH T ■ T ' ma9e S ' 9nal ' SU5 " band HLl inClUdSS the bi9h-frequency components a.ong the fi« 
direction and the low-frequency components along the second direction. Finally, the sub-band HH, includes the high- 
frequency components along the two directions 9 

fourth ba S nd C s 0 L n L Th^ ^2 ■ ^ ^ SUb " band ° f '° W ' rec < uencies LL i so a s, in the same way, to supply 
hen fiH^^h ^ ? 2 ,ntermediate level of resolution ^ the decomposition. A third analysis block 
then f.l ers the sub-band of low frequencies LL 2 so as to supply four sub-bands LL 3 , LH 3 , HU and HH a Finally in this 

HL 4 and HH 4 with the lowest resolution in this decomposition. 4 4 ' 

[0143] Hence thirteen sub-bands and four levels of resolution are obtained. Obviously, the number of levels of res- 
olution, and consequently of sub-bands, may be chosen differently 

K™L f >bViOU , Sly - ° ,her ,ypes of spec,ral transformation could be used, such as the discrete Fourier transform the 
discrete cosine transform or the Fourier-Mellin transform «"*iurm, me 

mli 5 L l be^nsertea. Way ' ' reqUenCy SUb " b3ndS M9 3 S6t ° f SpeC,ra ' coefficien ts into which a water- 

This sub Iand h SH f f Z enCy T^, HH ' " be conside ^ here in order to insert a watermark. 
This sub-band HH, thus supplies a set of coefficients X, equal to N in size, for example 

JJw - «! X = {X| ' 1 * ' S N} be the Set ° f coefficients representative of the digital image I 

S i. ^ * *° ' SUbSet ° f ° f thiS S6t X " representative on.y of a part 

Sa^dom JZTh l ° inSe H rt 3 T termark ,hiS Set ° f coefficten ts X, one technique consists in inserting a pseudo- 
5£T TK fw T< S l PeCUUm S ° 38 l ° fender thiS Si9nal unde tectable by spectra, or statistical analysis, 
averse of lerji2 m I ? 3 "T™* ° f SiZe P ~ K * P seud o-random signal with known distribution and 
11 S 1 T 7 des Pread distributions are the binary distribution {-1,1}, the uniform distribution over [-1 , 

1] and a centred normalised Gaussian distribution N (0. 1 ). 
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[0151] The linear insertion lormula is : 



X'j = Xj + b ttj w r with 1 < j < P, 

5 

in which {Xj, 1 < j < P} is a subset of the set of coefficients X, b is an information bit, and is a weighting coefficient, 
also called modulation amplitude. 

[0152] The set {j, 1 < j < P} is also called information-bit insertion carrier. 

[0153] In order to insert a signal which may be formed from several information bits, it is worth partitioning the set 
w of coefficients so as to determine the number of regions or the number of insertion carriers available in this set. 

[0154] To do that, the insertion device 10 in accordance with the invention comprises a device for determining a 
partition 11 and insertion means 12. 

[0155] The device for determining a partition 11 is suitable for carrying out adaptive partitioning of the set of coeffi- 
cients as a function of the detectabiiity of an inserted watermark. 
is [01 56] This method of adaptive partitioning consists in varying the number of coefficients modulated in order to insert 
an information bit of the watermark, so as simultaneously to satisfy the criteria which are set of invisibility and of correct 
detection probability. 

[01 57] A description will first of all be given, with reference to Figure 3, of a method of determining a partitioning into 
distinct regions of a digital image which makes it possible to maximise the number of regions obtained in order to 

20 maximise the number of information bits inserted in the image. 

[0158] Consideration is given to a representation of the digital image by a set of coefficients, either in the spatial 
domain, or in a transformed domain, the coefficients being in the latter case hybrid, that is to say located both in the 
spatial domain and in the frequency domain. Such a representation of the image is for example obtained by using a 
sub-band issuing from a spatio-f requency decomposition of the image, as a discrete wavelet decomposition previously 

25 described. 

[0159] The partitioning of the image is effected in order to insert a watermarking signal by modulation of the coeffi- 
cients of each region. 

[01 60] According to the invention, this partitioning into regions is performed by an adaptive partitioning as a function 
of a criterion of detectabiiity of an information bit inserted on this region. 
30 [0161] The number of regions obtained under constraint of the detectabiiity criterion makes it possible to define a 
maximum capacity of the digital image. This maximum capacity corresponds to the number of information bits which 
it is possible to insert in the digital image whilst complying with the detectabiiity criterion. 

[0162] Initially, at a step E100, this capacity Q is initialised to 0 and the starting point is an initial partitioning which 
is limited here to a single region R of size P. 
35 [0163] Naturally, the initial partitioning could be different and notably include already several distinct regions with 
sizes which might be different. 

[0164] It is assumed that the size P of this initial region R is greater than a minimum size denoted P stat which corre- 
sponds to the minimum size of a statistically significant sample of coefficients for detecting a watermarking signal 
inserted in these coefficients. Typically, this minimum size can be fixed at 100. 
40 [0165] This initial regbn R can consist of the set of coefficients of the representation under consideration or a subset 
of this representation. 

[0166] A test step E101 then tests the detectabiiity of the insertion of a watermarking signal in the initial region H. 
[0167] This detectabiiity test can be effected by actual modulation and detection operations on the coefficients of 

the initial region R. . 
45 [0168] In practice, all the coefficients of the initial region R are copied into a working memory so as not to directly 

watermark the image. 

[0169] Next all the coefficients of the initial region R are modulated by inserting the watermarking signal. 
[0170] Conventionally, the following procedure is followed in order to effect this modulation: 
[0171] The initial region R consists of a set of coefficients: 
50 X: {Xj, i= 1, ...P) 

[0172] Normally a linear modulation is used of the type. 

X i * = X l ±a i W j , withi = 1,...,P 



55 



where W : {W j: i = 1, ...,PJ is the watermarking signal generally chosen as a pseudo-random sequence with a 
known distribution and a nil mean. 

[0173] The value of the sign ± of the modulation depends on the binary value to be inserted: for example the sign - 
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corresponds to the value 0 and the sign + corresponds to the value 1 

thelnte Jalt-riT' 1316 ^ ChO0S6 ^ ^ " * pseud °- random sequence W which follows a uniform law on 

f°c 1 I?l f^T"* f ny P seud °- random watermarking with known distribution and nil mean can be suitable The most 
usual distributions for the watermarking signal W are, apart from the aforementioned uniform distribution [-11 the 
binary d.stnbut.onM,1> and the centred standardised Gaussian distribution N(0 1) ™*noui.on [ i, i], the 

[01 76] The modulation can optionally be protected by means of a secret key characterising for example the kernel 
for reproducing the pseudo-random sequence W. 

EiUS ThG tGrm % ' S the modula,ion amplitude, or weighting coefficient, applied to the coefficient X 
[0178] A constant weighting coefficient can be used for all the coefficients of the region R, such that a- = o v for anv 
., the value of the weighting coefficients «, being equal to a maximum value ensuring impercep.ibility of the waLmark 
ing signal for the modulation of the set of coefficients representing said digital image 

thll? th Th ;f h maximum va,ue "v < referred to a * ^st Noticeable Difference or JND amplitude) decreases in fact with 
m«m 9 r ^T"" watermarking signal W, and consequently with the number of modula.able coefficients 
al ZLTZ h © capacity Q will be maximised when all the coefficients of the representation of the digital image 

coeff^n ,T°! , * ' S )Ud ' CiOUS t0 dGtermine 8 modulatlon JND amplitude for a number of modula.able 

coefficients equal to the cardinal number of the set of coeff.cients representing the ,mage 

[0181] This maximum value a v or JND amplitude corresponds to the maximum modulation amplitude which can be 
used .n the aforementioned linear insertion model beyond which an observer is capable of visually detecting a change 
in the reconstructed watermarked image. y 9 

SIS J" ° rd KT t t0 d f e,ermine tWs JND am P'^e, it is possible to use a visibility model which makes it possible to 
Inti h visibility of a watermarking operation according to different parameters, namely the representation of the 
h To mh 9 3 Spa,,0 - fr f et ? uenc y transformation S used, the sub-band considered for the insertion the type of distri- 
bution of the sequence of the watermarking signal W and the length P of the sequence W 

noise- A B S WATSON 2 !nS1 * """T" * ^ ^ Sntit,ed " Visibilit y of wave,et quantization 

f WAT S° N et al - ^ Tr ans. on Image Process, 6(8), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of modulated coefficients from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for the detailed description of this model 

[0184] Thus use is made of a function dependent on the length P of the watermarking signal W, the type of trans- 
formation S used and the (base) sub-band considered for the insertion, but independent ofall the coefficients to be 
watermarked X. 

[0185] This calculation function used can be written: 

ol v (P > S, w) = 



base^ 



where 



<W(S) is a base value, dependent on the transformation used S and on the (base) sub-band considered for the 
insertion of the maximum weighting coefficient ensuring imperceptibility during modulation of a single coefficient 
or tnis sub-band, 
|3 is strictly greater than 2, and 

E[MP] is the mathematical expectation of the function IwlP. 

h°o!^ 61 T v 6 b3 f 6 V3lUeS °' base(S) ° an be measure d once and for all, for each sub-band of coefficients in the wavelet 

™ ^ PC T ' m 3 S ' ngle P s V chovisual measurement, and be stored in a table of visual amplitudes 

[01 87] p is the exponent of a Minkowsky sum and can be chosen equal to 5 for example 

H ™ e ^.ematical expectation £[I^P] corresponds to an estimation of the mean of the function MP 

Th i!fi T S \ . y m0d , e ' d ° eS t3ke im ° aCC ° Unt the image 1 to be watermarked itself and is independent thereof 

Th s ,s equ,valent to considering that the image . is uniform. It is a question of a "worst case" model since the presence 
of the image signal makes it possible to visually mask the modulation itself 

[01 90] The JND amplitude can thus be calculated for a length P of the watermarking signal W equal to the cardinal 
number of all the modulatable coefficients of the sub-band under consideration qua. to tne cardinal 

TK° WeVer ',i! ^ ' UdiCi T '° 6Xpl0it thS <aCt th3t ,he image Si9nal itself makes « P° ssible *> <he modulation 
i™L ^ ■ e ?T P ' USS iS mad6 ' f ° r eaCh coefficien < *i <° be modulated, of a weighting coefficient in 
accordance with a law of the form c, = I, c^, where k f is a modulation factor dependent on the coefficients situated 
close to the coefficient considered on the region R and a,, is equal to the maximum JND value of the weighting coef- 
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ficient. 

[0193] Thus each coefficient will be modulated according to the local content, which makes it possible to locally 
increase the amplitude of the modulation to the benefit of detectability. 

[0194] In order to test detectability in the test step E101, a detectability amplitude is next calculated from the de- 
5 tectability criterion and this detectability criterion is validated by comparing the detectability amplitude with a predeter- 
mined threshold value denoted Tc. 

[0195] A calculation of correlation C(X*, W) between all the modulated coefficients X* and the marking signal W is 
generally made and it is decided whether there has actually been watermarking of the image when the result of the 
correlation calculation C(X* : W) is greater than the predetermined threshold value Tc. 
io [0196] A standardised statistical test can also be used for the detection, such as the one described in the article 
entitled "A method for signature casting on digital images" by I. PITAS, in Proc. ICIP, pages 215-218, September 1996. 
The detection is then characterised in terms of probability. It is thus possible to choose a threshold value Tc corre- 
sponding to a fixed detection probability level, for example 99.95%. 

[0197] Thus, for a given region R defining a set of coefficients X of size P and with a fixed weighting coefficient a, it 
15 is possible to calculate a detectability amplitude and compare it with a.threshold value Tc so as to validate or not the 
detectability criterion on this region R. 

[0198] The detectability criterion is in fact in theory a function depending on the size P of the region R, the detection 
threshold Tc related to the required detection probability, the watermarking signal W, the variance a 2 x of the coefficient 
X to be modulated on the region R and the weighting coefficient used a for modulating the coefficients X according to 
20 the linear model described above. 

[0199] Thus, for example, for a detection by standardised hypothesis test and a value of the weighting coefficient a 
which is constant and equal to a v for all the coefficients X { of the region R, a minimum length P min (X) is defined by the 
function; 

Pmin(>0 = (^f) [t>o 2 x + co 2 + d] 

where a, b, c and d can be assimilated to constants The signal being non-stationary, that is to say the variance 
30 a 2 x being a function of space, it is indeed necessary to adapt the modulation length P generated by the watermarking 
signal to the image signal itself. 

[0200] On decoding, the set of coefficients X being modulated, this function also depends on the variance a 2 w of 
the inserted watermarking signal W. 

[0201] In addition, it is necessary to provide for any distortions which the image may undergo after modulation, due 
35 for example to a compression of the image. 

[0202] When this noise can be considered to be additive and decorrelated from the coefficients, it is possible to use 
a model dependent on the variance o 2 q of this additive noise modelling any distortions applied to said coefficients after 
modulation. 

[0203] Thus, to a given weighting coefficient a, making it possible to ensure invisibility of the inserted watermarking 
40 signal, a detection threshold Td fixed as a function of the detection probability required at decoding and a set of selected 
coefficients X, there also corresponds a minimum modulation length P min (X) which makes it possible to ensure de- 
tectability of the inserted watermarking signal. 

[0204] On decoding, for a detection by standardised hypothesis test, this minimum length P mjn (X) is then defined by 
the function: 

45 

(9Tci\2 2 2 2 2 2 

Pmin(X)=[^J [6(a x +cc o £ w +c q ) + ca + cf] 

so where a, b, c and d can be assimilated to constants. 

[0205] It is therefore necessary, if the same detectability test is used during coding and decoding, to choose a thresh- 
old value Tc on coding greater than that Td used at the time of decoding if it is wished to be able to find the same 
partitioning on decoding of the image. 

[0206] In practice, the threshold value Td on decoding which it is wished to use is fixed, that is to say for example 
ss the required detection probability level when the detection is made by a standardised hypothesis test, and a higher 
threshold Tc on coding is chosen, in order to take account of the variance in the watermarking signal W and the variance 
of an additive noise fixed a priori so as to anticipate the influence of the subsequent processing undergone by the image. 
[0207] This condition thus makes it possible to find the same detectability criterion which constrains the partitioning 
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of the digital image without having to transmit to the decoder the partitioning effected at the time of insertion of the 
watermarking signal W. 

[0208] In this example embodiment, a partitioning by iterative division of the digital image is effected It is conse- 
quently considered, at step E101, that the detectability criterion is satisfied when the calculated detection amplitude i* 
greater than or equal to the predetermined threshold Tc. 

[0209] If the detectability criterion is not satisfied, then it is probable that the detectability criterion is also not satisfied 
in subregions of the initial region R of smaller size because the detection probability generally decreases with the size 
of the regions. A step E102 of stopping the method of determining a partitioning makes it possible to stop the process 
of iterative division of the image. 

[0210] Insertion not being judged possible in the initial region R, the capacity Q of the image is therefore nil 
[0211] If on the other hand the detectability criterion is satisfied, that is to say here that the detection amplitude i- 
greater than or equal to the threshold value Tc, it is possible that it is also satisfied for subregions of the initial region R 
[0212] It is therefore judicious to attempt to divide the region R as long as the detectability criterion is satisfied in 
order to increase the capacity Q. 

[0213] There are of course several possible means for effecting the spatial division of the region R into subreqion* 
Rj, of sizes Pj, with j = 1 , ... , K. 

[0214] In this example a division into a quadtree is chosen, as illustrated in Figure 4. 
[0215] This division consists of dividing each region into four subregions of identical size. 

[0216] Thus, by way of example, in Figure 2, a digital image of initial size N is being partitioned into two sub-region- 
R5 and R6 of sizes N/4, seven subregions R1, R2, R3, R4, R7, R8, R9 of sizes N/16 and four subreqions R10 Rn 
R 1 2, R1 3 of sizes N/64. 

[0217] Thus, at the division step E103, the region R is divided into subregions R, j = 1,.... K of sizes P, In a case of 
a division into a quadtree, K = 4 and Pj = P/4. J 
[0218] At step E104 a first subregion of index j (for example j = 1) is considered and a value of potential gain in 
capacity denoted G Q is initialised to - 1 . 

[0219] The partitioning determination method then includes a prior step E105 of comparing the size P, of the region 
Rj with the minimum size P stat corresponding to the minimum size of a statistically significant sample for the detection 
of a watermarking signal inserted in said region as described above. 

[0220] If the size Pj of the region R. is less than this minimum size P stat , at step E 1 08 the index J = J+1 is incremented 
in order to process the following sub-region Rj +1 . 

[0221] Otherwise, the detectability criterion is checked on the subregion R= in the same way as previously described 
with reference to step E 101 . 

[0222] In practice, the watermarking signal is inserted on the coefficients of the said region Fl, a detection amplitude 
is calculated and this detectability amplitude is compared with the threshold value Tc. 1 
[0223] The detection amplitude can be a correlation value between the modulation signal X* and the carrier of the 
watermarking signal W or the result of a standardised hypothesis test. 

[0224] If the detectability criterion is verified, the value of potential gain in capacity G Q is increased by one unit at 
step E 107. 

[0225] Then at step E108 the index j of the subregions is incremented and it is checked at the test step E109 that 
this index j is less than or equal to the number of K. 

[0226] In the affirmative, steps E105 to E109 are reiterated on the following subregion. 

[0227] Otherwise, when the index j = K+1 , it is checked in a test step E11 0, whether the potential gain in capacity 
value Gq is strictly positive, that is to say whether there exist at least two subregions R of the initial region R for which 
the detectability criterion is validated 

[0228] I n the negative, the partitioning of higher rank carried out on the initial region R is not validated and this region 
R is preserved. Such is the case when the detectability criterion has not been validated for any of the subregions Ft 
or only for a single subregion Rj. ' 
[0229] In the latter case, it is preferable to keep the region R for inserting the watermarking signal since it is of size 
P greater than the size of the single validated subregion Rj. 

[0230] At step E112, the value of the capacity Q is then increased by one unit and the region R of size P in the 
partitioning is kept. 

[0231] An actual insertion step E113 by modulation of the coefficients of the region R can then be implemented in a 
normal fashion using the linear insertion model described above. 

[0232] An elimination step E114 then makes it possible to eliminate the region R from a stack of regions to be proc- 
essed. 

[0233] On the other hand, if at step E110 the value of the potential gain capacity Gc is strictly positive, in an addition 
step E1 11 all the subregions Rj for which the detectability criterion is validated in the stack of regions R to be processed 
are added. 
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[0234] A return step E115 makes it possible to reiterate all the steps E103 to E11 3 in order to process the regions 
of the stack in a recursive tashion. 

[0235] At the end of the recursive partitioning method just described, a capacity value Qm corresponding to the 
maximum capacity of the image for insertion of a watermarking signal is obtained, guaranteeing a correct detection of 
this signal. This maximum capacity Qm is equal to the cardinal number of the distinct regions obtained by the said 
partitioning on which the detectability criterion is validated. 

[0236] In addition, the insertion step proper E113 can be omitted during a first "pass" which then makes it possible 
only to determine the maximum capacity Qm of the image and the validated regions. 

[0237] Once the maximum capacity Qm has been calculated and the corresponding partitioning is stored in memory, 
it is then possible for a user to determine the message which he wishes to insert, corresponding to a number of infor- 
mation bits equal to the maximum capacity Qm, and then to insert this message in a normal manner by modulating 
the coefficients of the different regions obtained by the previously determined partition. 

[0238] A description will now be given, with reference to Figure 5, of a second embodiment of the invention in which 
it is wished to insert a message of fixed capacity Qs in an image I. Naturally, this fixed capacity Qs must have a value 
less than the maximum capacity Qm of the digital image to be watermarked. _ 
[0239] As illustrated in Figure 5, the method of determining a partitioning in this embodiment includes steps E100 to 
E113 identical to those of the method of determining a partitioning of the embodiment of Figure 3 and which it is not 
necessary to describe again here. 

[0240] As before, at step E 11 2, the value of the capacity Q is incremented when the partitioning of rank directly higher 

is not validated. . . ... 

[0241] Then the method includes a comparison step E11 6 in which this value of the capacity Q is compared with the 

fixed capacity value Qs. 

- [0242] In practice, it is tested at step E116 whether the capacity value Q is equal to the fixed capacity value Qs. 
[0243] In the affirmative, a stop step E117 makes it possible to stop the iterative partitioning of the image since there 
are a number of regions equal to the number of information bits Qs which it is wished to insert. 
[0244] Otherwise, in an elimination step E114, the processed region is eliminated from the stack of regions R to be 
processed formed at step E11 1. 

[0245] The method of determining a partitioning also includes, in this embodiment, a sorting step E118 which sorts 
the regions R to be processed in the stack by detectability amplitude value calculated on these regions R. 
[0246] Thus a return step E115 reiterates the partitioning from the division step E103 for a region R in the stack 
having the highest detectability amplitude amongst all the other regions to be processed. 

[0247] A sorting variant can consist of sorting the regions to be processed in decreasing order of size so as to process, 
as a priority the largest regions which a priori provide better detectability. 

[0248] Naturally, a partitioning by iterative merging of the digital image could be carried out in the same way, for 
example also using a quadtree partitioning structure. 

[0249] The initial partitioning then consists of regions with smaller acceptable sizes, greater than or equal to the 
previously defined minimum size P stat . 

[0250] For each region of the partitioning for which the delegability criterion is not validated, this region is merged 
with another region of the partitioning. . 
[0251] This other region of the partitioning is if possible a region for which the detectability criterion is not validated 

so as to increase the capacity. ^„„ ti 

[0252] Alternatively, another region can be a region satisfying the detectability criterion so as to make detection 

rolidblo 

[0253] ' There are many possibilities for the choice of the region to be merged. It is possible to take as the merger 
criterion the detection amplitude calculated on each of the regions and merge for example the regions associated with 
the lowest detection amplitudes. . . „ ^^o.^h 

[0254] In this particular case of a merging in a quadtree partitioning structure, the merging is structurally constrained 
and consists of deciding whether or not tour adjacent blocks, for example R1, R2, R3 and R4 in Figure 4, must be 

[OMsf A description will now be given, with reference to Figure 6, of the associated decoding method which makes 

it possible to find the partitioning when the watermarking signal inserted in a digital image is extracted. 

[0256] Indeed, because of the adaptation of the partitioning to the digital image, this partitioning will not be known a 

priori at the decoder. . 

[0257] The decoding method then consists, in general terms, of testing all the acceptable partition mgs. The same 
iterative implementation technique is used as the one used when the partitioning is determined for inserting the wa- 
termarking signal. ... . 
[0258] Thus a description will be given of an example embodiment in the case of a segmentation by dividing by 
blocks in a quadtree, symmetrical to that using the partitioning determination method of the embodiment described 
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with reference to Figure 3. 

S. „1' S " P E20 °' ,s ,n """ P»n»^'"9 Med lo a S in flte ,. 9ksn R o, p TWs ,.„,„ 

[0262] In practice, a detection amplitude is calculated for all the coefficients of the region R and the detectabi.i.v 
SST '.t rooSe ^ C ° mpar ' n9 h thiS de,eC,ion am P |ltude with a predetermined threshold value Id or decS ng " 

EdS2S eed by means of a s,andardised statistical ,est ' and ,he ™* c -p- d ^ *» 

th ^ 9 ° f ,he ' ma9S 6Ven When the la,ter has been mad * no^y by post-processing 

hat .to say, ,n th.s example, a division of the region R into four subregions R, of identical sizes. 

Sm sleo E2u 2 T e HrH SteP SUbregi ° nS * ^ added h 3 S,aCk ° f re 9 ions R to be P™~— and, a, 

a return step E202, the decodmg method » rerterated recursively as from the test step E201 on each region R of the 

S22 JUT 31 I* 16 . StSP ° f Ca J CUlatin9 ,he det ectability criterion E203, the latter is verified for a region R there is 

fn Z dio TT S,6P E20S information blt * this region R as an element of the method inserted 

esu o Hn J tTst ^n rT!' 0n 1 e " e T " 3 COnVentional < ash ™ b V -Iculating a detection ampmude ^he 
result of the test gives, in absolute value, a detection amplitude to be compared with the threshold value Tri of tho 
decoding, the sign giving the value 0 or 1 of the information bit inserted 

ST? l pon°o Uld ^ n ° led ,hat thiS extrac,ion is a,rea dy partly effected during the detection step E203 so that the 
detection E203 and extraction E206 steps are structurally coupled 

[0270] An elimination step E207 is implemented in order to eliminate the region R from the stack of reaions to he 

X^SnZ C ^7 f ° rmed 1 addi, " n St6P E2 ° 5 ' and ,hS retUm Steb E202 toen makes t possibt to ^tlnt 
the decoding method recursively on each region R of the stack to be processed implement 

[0271] A description will now be given, with reference to Figures 7 and 8 of a practical aeolication of the moth ^ of 
^ermining a partitjonjng of a djgita| jmage jn g pgrtjcu|ar em 9 bodjment apP " Ca, '° n of ,he method °< 

o S app ' ication ' the watermarking signal is inserted by means of a technique of insertion by spectrum 
oulnoJt 9 ' f m °T at, ° n ? co ^«=ients of a space-frequency representation of the image, obtained by a spSre 
q^ency transform*™ S of the digital image I as prev.ously described with reference to figure 2 P 
[0273] In genera, terms Sj sub-bands are obtained, with J = 1 . .... M and M equal here to 13 at a transformation step 

\ZlentZ\^^T^T S) ° f indSX j C ° nSidered 31 St6P E3 ° 1 • 3 partitionin 9 is ef,ec,ed a < step E302 by 
mplementing the method of determining a partitioning described previously to with reference to Figure 3 

[027S] In this way a capacity Qj is derived for each sub-band S- 

^ regtn^vrdalr^ 0 " * m ° dU,a, '° n °' C ° eff ' Cien,Js ™ makeS » P ° SSib,e to '" sert Q i ^ion bits 
[0277] m steps E304 and 305, the following sub-bands are next considered as long as all the sub-bands S have not 

szK^srs;. E302 and E3 ° 3 ot partitionins and — 9 a -~ -rsas^s 

Stv SeouanLT 56 ,fan f t ° h rmati0n S ° is a PP' ied to toe image in order to obtain the watermarked image, of total 
r^a V q SUm ° f the ca P ac,ties Q , calculated for each of the sub-bands 8, 

s ou ?^cXtn n d d :;^ 

EHl^c! , !;^l^ w ,h 9 T^ 9 T th0d " ,hiS embodi — is implemented in a similar fashion using 

srcri:r ce 9ure 6 ' the initiai re9k>n r here correspondin9 each time to a 

K J fTr!f ^c' 03 ' '° S,6P E3 °°' diVid6S ima9S ' int ° S i SUbband * With J = 1 .... M 

[0282] A first sub-band Sj ,s considered at step E401 The order of processing of the sub-bands is identical to that 
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used during coding in order to extract the inserted message in order. 

[0283] The adaptive partitioning is implemented at the partitioning step E402 in order to calculate the capacity Q, of 
the sub-band Sj. Next, at the extraction step E403. the Qj information bits are extracted by demodulating the coeff.c.ents 

of the Q regions of the partitioning 

5 f0284] At the test steps E404 and E405. it is checked in the usual fashion that all the sub-bands have been processed 
and. in the negative, the decoding steps E402 and E403 (partitioning and extraction) are reiterated on the remam.ng 

[0285] nd |n this way. at the output of the decoder, a message of Q bits is obtained, Q being equal to the sum of the 

capacities Q: of each sub-band Sj. 

10 [0286] In a variant of this method illustrated in Figure 9, it is possible to effect a partitioning on a post-processed 
image signal in order to obtain an a priori robustness of the watermarking signal to any post-process.ng undergone by 

[0287] a9 The principle of this variant consists of performing a detectability test in a test step E309 which is no longer 
coupled to the steps of partitioning E302 and insertion E303 described previously. „„ orf ™ m0 H 
is [0288] After the insertion of the watermarking signal in a region, a reverse trans.ormat.on S-i .s first of all performed 
at step E306 in order to find the watermarked image. Then, in a step E307, post-processing is applied, such as a 
compression of the digital image. . . h , 

[02891 Next in a step E308. the transformation S is reiterated in order to find the watermarked coeff.c.ents. possibly 
made noisy by the post-processing carried out on the image at step E307, and the detectability criterion is calculated 

20 as previously described. . „_„„i^ k., 

T02901 Steps E300 to E309 are reiterated on another acceptable partitioning of the .mage, here for example by 
dividing the relevant region R into four, and the detectability test is performed at step E309 on the four sub-regions. 
[0291] This partitioning is adopted only if it represents a gain in capacity, that is to say. in th.s example of quadtree 
division, if the detectability criterion is validated on at least two sub-regions. 

ss [0292] This variant embodiment makes it possible to directly take into account the post-process.ngs wh.ch w.ll be 

aoDlied to the diqital image. .., , . 

[0293] A method of determining a partition according to a preferred embodiment of the invent.on w.ll now be de- 
scribed, by reference to Figure 10. 

r02941 The method of inserting a watermark is employed on a digital image I. 
30 [0295] This insertion method first of all includes a method of determining a partition of a set of coeffic.ents represent- 

[0296] f 1 Thfs' partition^applied as a function of a criterion of detectability of an information bit inserted on each region 

raWTTlt should be noted that, in practice, a partition of a digital image into a maximum number of regions satisfying 
3 s a detectability criterion can be sought so as to ascertain the practical capacity of this image with regard to he insertion 
of a watermark. This maximum capacity of the digital image corresponds to the number of information bits wh.ch it « 
possible to insert into the digital image while complying with the detectability criterion. 

[0298] It is also possible, for a given watermark, comprising a given number of information b.ts of less than the 
maximum capacity of the image I, to seek the best possible partition of the image I which makes . poss.ble to insert 
this watermark of given length, for a given detection probability, and in regions of the largest poss.ble sizes so as best 

to ensure invisibility of the inserted watermark. ce nn-nh.«i««i.. 
[0299] The method of determining a partition first of all includes a step of spectral decomposition E500 of the wavelet- 
decomposition type described previously (Discrete Wavelet Transform or DWT). 

[0300] Thus a set of spectral coefficients X is obtained, representative of the digital image I. These spectra coeffi- 
45 dents are distributed into frequency sub-bands as illustrated diagrammatically in Figure 2, on completion of the de- 
composition step E500. The size of the set of coefficients X is equal to N = M x M for a square image. Obviously, the 
same method could be applied to a rectangular image. rQ „ H ^ m 
[0301] In this embodiment of the invention, the method includes a step E501 of generating a centred pseudo-random 
sequence of size N equal to the cardinal number N of the set of coefficients X representative of the .mage I. 
so [0302] Let w = fw mn , 0 <, m. n < M} be this pseudo-random sequence. 

[03031 This pseudo-random sequence w is formed from centred pseudo-random sub-sequences w*. 
[0304] The pseudo-random sequence can be seen as the merging of the pseudo-random sub-sequences w. 
[030S] Here, the centred pseudo-random sequence w is formed from pseudo-random sub-sequences w* of ident.cal 
size which correspond to a block, called unit block, of the set of coefficients X. 
55 [03061 Obviously, the pseudo-random sub-sequences w* could differ in size from one another. 

03071 Here the size of these unit blocks corresponds to the minimum number of coefficients suitable for formmg a 
reqion of a size which is statistically significant for satisfying the detectability criterion over th.s region. 
[0308] As described later, when a standardised statistical detection test is applied to a set of coefficients, a minimum 
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number of coefficients is necessary for this test to be significant 

number ^SS^T^ + -» ^ ^ than or equa, to this minimum 

[0310] For preference, pseudo-random sub-sequences w* are chosen to be equal in size to this minimum size This 

fror a ?inS, ' S ^""r'f ^r^ 5 Wh8n Partl,i ° n iS aChi6Ved * ««*^ combina?^ 
rom an initial segmentat.on of the image into unit blocks of size L x L The initial segmentation then cor espond -to 

onger statistically significant for the insertion and detection of an information bit generated are no 

[0311] For example, L = 8 is chosen. 

K PSeUd °" rand ° m SequenCe w is then the ™VW of N/64 pseudo-random sub-sequences w* forming unit 

[0313] On each unit block, in accordance with the invention, a centred pseudo-random signal of predefined distribu- 
tion (Gaussian, un.form, etc.) has to be generated on the basis of a function depending on a secret kev K and or ■ «"» 

of Te si°nn e a ? C wl ■ ^ ™ S ^ ^ * f (K ' k) ' « is ,hen necessa ^ * —-that n P "aJS ^he average 

of the signal w* is exactly equal to zero over each unit block average 

[0314] A first method for generating such a centred pseudo-random sequence consists in centring the sequence in 
a determ.n.st.c way, by symmetrica This method is suitable for generating a pseudo-random .uSiJS^ 

[0315] in this case, each centred pseudo-random sub-sequence w* is created by generating one half A/2 of the 
pseudo-random samples of the sub-sequence via a pseudo-random number generator of known centred d^ibu ion 
lle V s7 b ~! wl 6 P ~ d ° m «1 aerated in order toLain the other ha^t^es 

[0316] For a single-dimensional sequence, a set of samples of size A/2 { Wj , 1 < j < A/2} is therefore generated and 
his set ,s made symmetric in order to obtain the other half of the samples d = -w A . A/2 < j < A) 9enera,ed a " d 

E J2 Jl ° fd T l ° r h e t in, °T, ,he rand ° m Ch8raCter °' the Sec < uence thus obt ^ ed - Permutation can be carried out on 
the set of samples obtained { Wj! 1 < j < A} by using a secret key 

EndSn J? r C ° nd ? 96neratin9 a centred Pseudo-random sequence consists in creating each centred pseudo- 

S2^^2 UenM 9enSra,in9 PSeud °- random sam P les vb a pseudo-random number generator and by 

m^I ° f ,h6Se P seud °- ra ndom samples generated over each sample Y 

S th ^ J' 81 0, , al1 ' ^ USin9 S6Cret k6y f(K ' k) ' 3 s *W e ™e °< pseudo-random samples of the same size A 
as the unit block, i.e. {w h 1 < i < A}, is generated. 

[0320] Next the exact sum of these samples is calculated: 




[0322] Thus wk = (Wj-Sp/A, 1 S i < A} is obtained 

a°wa!\haTth7Z^n and T "TT^ £ iS ,hUS C ° mP ° Sed ° f ,he me ' 9in9 of the centred sub-sequences w* in such 

Lb 3 !,. « P art,,, ° n,n 9 of a se < °f coefficients modulated by the centred pseudo-random sequence w gives rise to 

subsets of coefficients also modulated by centred pseudo-random sub-sequences 

£22 ° bV Tf ly ' ° ,her ,echni q"es could be used to generate centred pseudo-random sequences 

ieauence w ^ ' ?' * E502 °' mOdU ' ati0n ° f the Set 0f ^efficients by the centred pseudo-random 

rzr;r men,8d so as to insert an information b » ° nt ° ,he - * s 
i 0 cZuti~r o rt ssss: ssr into a workina memo ^ ,or exampie a random - access — - 

KwJ!^ W ^ ,0n iS m ° dUlated ° nt ° S6t ° f C ° effiCientS ' f ° r —P' e b - 1 according to the 

X'j = Xj + a s Wj, with 0 < j < N. 
[0328] The weighting amplitude a, is chosen in the usual way so as to guarantee invisibility of the inserted information 
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[0329] A constant weighting coefficient can be used for all the coefficients, such that a. t = a v for all i, the value of the 
weighting coefficient c* v being equal to a maximum value ensuring imperceptibility of the watermark for modulation of 
the set of coefficients representative of the digital image. 

[0330] It may be worthwhile exploiting the fact that the image signal itself makes it possible to mask the modulation. 
5 [0331] To this end. it is possible., for each coefficien* Xj to be modulated, to use a weighting coefficient according to 
a rule of the form a t = kj.o^, where k t is a modulation factor depending on the coefficients situated in the vicinity of the 
coefficient in question and a v is equal to the maximum value of the weighting coefficient ensuring imperceptibility of 
the watermark for the modulation of the set of coefficients. 

[0332] Thus, each coefficient will be modulated as a function of the local content of the image, which makes it possible 
io locally to augment the amplitude of the modulation to the benefit of the detectability. 

[0333] For preference, in this embodiment of the invention, it is wished to apply a partitioning to a post-processed 
image signal so as to obtain robustness, a priori of the inserted watermark, against any post-processing undergone 
by the image I. 

[0334] To do that, a distortion is applied, in a distortion step E503, to the set of modulated coefficients. 
is [0335] By virtue of the modulation of the coefficients in a single operation, on the basis of a centred pseudo-random 
sequence, the distortion step E503 can be implemented once only upon determination of a partition, in contrast to the 
prior techniques in which each sub-region of coefficients was modulated independently of the other and required dis- 
tortion to be applied to each level of the partition. 

[0336] This is all the more advantageous since, here, the step E503 of applying a distortion includes the following 
20 sub-steps as illustrated in Figure 10: 

inverse spectral transformation E503a of the set of coefficients X'j modulated in order to obtain a watermarked 
image; 

application proper, E503b, of distortion to the watermarked image; and 
25 - spectral transformation E503c of the watermarked image in order to obtain a set of modulated coefficients after 
distortion. 

[0337] The distortions applied are distortions which alter the value of the pixels without altering the geometry of the 
image, such as the addition of noise, alteration of the contrast, compression, low-pass or high-pass filtering, for exam- 
30 pie. . 

[0338] Compression of the image is considered here, according to the JPEG standard (for Joint Photographic Expert 
Group) with specified quality factor with q = 75 typically by default for JPEG compression. 

[0339] In practice, an inverse wavelet transformation DWT" 1 is applied here, then compression followed by decom- 
pression of the image, and finally forward wavelet transformation again in order to be able to estimate the detectability 
35 _ of the information bit inserted in the spectral domain. 

[0340] Next an iterative partitioning of the set of coefficients is applied, so as to search for a partition of the set of 
coefficients in which each inserted information bit can be detected correctly according to a fixed detection probability. 
The number of regions thus determined during the partition corresponds to the practical insertion capacity of the image 
I. 

40 [0341] The means for displaying to the users Q bit determined by calculating the capacity should exist. 

[0342] Various techniques for iterative segmentation by division/combination can be used as previously described. 
[0343] A first way of performing the preferred embodiment of the invention will firstly be described, by reference to 
Figure 11, in which the partitioning is achieved by iterative division on the basis of an initial segmentation of the set of 
coefficients into regions of a statistically significant size for validating a detectability criterion over these regions. 

45 [0344] Following on from the steps ol transformation E500, generation E501 , modulation E502 and distortion E503 
as described above, a first frequency sub-band will be considered, in a step E504, as a set of coefficients to be parti- 
tioned. 

[0345] The initial segmentation of this set of coefficients corresponds to a spatial segmentation of the image I into 
blocks of unit size as defined above. 
50 [0346] By its very construction, the pseudo-random modulation sequence is centred on each of these unit blocks. 
[0347] In a quadtree-mode combination, it is always four blocks of equal size which are grouped together, called 
"child blocks" in order to form a "parent block", and this combination is recommenced iteratively, giving the notion of 
recursiveness. 

[0348] Obviously, partition takes place frequency sub-band by frequency sub-band of the signal arising from the 
55 spectral decomposition. 

[0349] Each unit block is considered here, in a step E505, in an order of scanning of the sub-band, for example from 
left to right and from bottom to top. 

[0350] A validity criterion is applied over a first unit block which has not yet been processed. 
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[0351] In a calculation step E506, a statistical detection test is applied to this unit block of modulated coefficients XV 
[0352] An example of a standardised statistical test is given below. 

[0353] Let X'j, 1 < i < P be the set of coefficients to be tested after insertion of the pseudo-random sequence w and 
application of pre-distortion. 

[0354] The test consists in calculating, for this set of size P: 



10 



15 



t=MJp 
Jv 

in which M represents an average value of correlation between the values of each coefficient X'j and of modulation 
Wj, and V the variance of these values. 
[0355] Hence in practice: 

M = - Iw^^-tfJ .with M x = 

p i-i ' p 
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[0356] The calculation of this test value T requires knowledge of the pseudo-random modulation sequence w which 
is generated once and for all at the above-described generation step E501 
30 [0357] The detectability is given by the comparison of this test value T with a certain threshold P p , which corresponds 
theoretically to a probability of correct detection of p% if the test T follows a Gaussian law with variance equal to 1 for 
the two hypotheses tested (presence or absence of watermark). 

[0358] A comparison step E507 makes it possible to compare this test value T with the threshold value P p or correct 
detection rate. 

35 [0359] On completion of this comparison, if the calculated value T is higher than the correct detection rate, the block 
tested is held in memory, in a memory-storage step E509, as being a region available for inserting an information bit. 
[0360] In practice, a capacity-calculating step E510 can be implemented in such a way that, for every block of coef- 
ficients tested which is suitable for having an information bit inserted, the value of the capacity is increased by one unit. 
[0361] In contrast, if the calculated value T is less than the correct detection rate P a block of higher level is selected, 

40 in a combining step E508, that is to say, here : a parent block corresponding to the combination of four unit blocks, and 
the steps E506 of calculating the test value T and of comparison E507 with the correct detection rate P p are reiterated 
over this new block. 

[0362] During this combining of child blocks in order to form a parent block, those child blocks are preferably chosen 
tor which the detectability criterion is not satisfied. 
45 [0363] For application of this detectability criterion, any standardised statistical test can be used for the detection, 
such as that described in the article entitled "A method for signature casting on digital images" by I. PITAS, in Proc. 
ICIP, pages 215-218, September 1996, which makes it possible to characterise the detection in terms of probability. 
[0364] A threshold value P p can be chosen, corresponding to a fixed level of detection probability, 99.95% for ex- 
ample. 

so [0365] Hence, for a given region defining a set of coefficients X of size P and with a fixed weighting coefficient a, it 
is possible to calculate a detectability amplitude and to compare the latter with a threshold value P p so as to validate 
or not validate the detectability criterion over this region. 

[0366] In a general way, the detection probability decreases with the size of the regions tested, such that, by com- 
bining regions in order to increase the size of the block of coefficients used for inserting a watermark, the value of the 
ss test T is increased over this new block. 

[0367] Hence, on completion of these steps of adaptive partitioning E504 to E508, an image insertion capacity can 
be determined, equal to Q bits, as can an appropriate partition of the image corresponding to this capacity Q 
[0368] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit. into the 
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various regions of the partition. 

[0369] Obviously, the step E520 of insertion proper of the watermark can be coupled with the determination of the 
partition, in such a way that each region available on conclusion of the detectability test E506. E507, can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients of this region. 
5 [0370] The insertion step E520 makes it possible actually to insert an information bit of the watermark, by choosing 
a coefficient value b = 1 or = -1 in the above -described modulation formula, the alternating of the bits conveying 
the content of the message inserted. 

[0371] A centred pseudo-random sequence of known distribution can be used, different for each region of the par- 
tition, and initialised by a secret key K\ as illustrated in Figure 1 
10 [0372] It should be noted that this secret key K' may, if appropriate, be identical to the secret key K used during the 
generation of the pseudo-random sequence w for determining the partition. 

[0373] Moreover, a centred pseudo-random sequence w identical to that generated at the generation step E501 
could be generated. 

[0374] A partition of an image I is shown diagrammatically by way of example in Figure 10, on conclusion of steps 
is E504 to E508, in which, for example, 25 separate regions are obtained making it possible to insert a 25-bit watermark. 
[0375] Obviously, any other partitioning technique could be used. 

[0376] A second way of performing the preferred embodiment of the invention will now be described, by reference 
to Figure 12, in which the partitioning is achieved by iterative division of the set of coefficients. 
[0377] On conclusion of the steps of transformation E500, generation E501 , modulation E502 and distortion E503 
20 as described above, a first frequency sub-band will be considered, in a step E511 , as set of coefficients to be partitioned. 
[0378] A quadtree-mode division will be considered here, in which one block, called "parent block" is always divided 
into four blocks of equal size, called "child blocks". 

[0379] Obviously, partition takes place frequency sub-band by frequency sub-band of signal originating from the 
spectral decomposition. 

25 [0380] Hence, in a step E51 2, an initial block corresponding to a first frequency sub-band is considered. 
[0381] This initial block is stored in memory in a list L of blocks to be processed. 

[0382] Next, in a step E513, a first block B to be processed, coming from the list of blocks to be processed L, is 
considered. 

[0383] A detectability criterion is applied on this block B in a checking step E514 
30 [0384] In practice, as previously described, the value T of a statistical detection test on this block of modulated 
coefficients is calculated in a calculation step E514a. This step is analogous to the calculation step E506 described by 
reference to Figure 11. 

[0385] The detectability is then given by the comparison of this test value T with a certain threshold P p , which theo- 
retically corresponds to a correct detection probability of p% if the test T follows a Gaussian law with variance equal 
35 to 1 for the two hypotheses tested (presence or absence of watermark). 

[0386] A comparison step E51 4b makes it possible to compare this test value T with the threshold value P p or correct 
detection rate. This comparison step E51 4b is analogous to the comparison step E507 described by reference to Figure 
11. 

[0387] On completion of this comparison, if the calculated value T is less than the correct detection rate P p , a check 
40 is made, in a test step E515, as to whether there are blocks remaining to be processed in the list L 
[0388] If not, the recursive partitioning algorithm is terminated. 

[0389] Otherwise the following block B of the list L is considered at step E513, and the checking step E514 is reit- 
erated. 

[0390] It, on completion of the comparison step E514b, the calculated value of the test T is higher than the correct 
45 detection rate P p , detectability is tested over each child block of this parent block B in a checking step E51 6. 

[0391] In practice, the parent block B is divided into four child blocks and the detectability criterion is checked on 

each of these blocks according to a step analogous to the previously described checking step E514. 

[0392] A lest E51 7 makes it possible to verify whether the detectability criterion is valid over at least one child block. 

[0393] If yes, a division step E518 is implemented in order to actually divide the block B into four child blocks, and 
50 an addition step E519 allows these child blocks to be added to the list of blocks to be processed L. 

[0394] The test and checking steps E513 to E519 are then reiterated over each of the blocks of the list L until the 

latter is empty. 

[0395] In contrast, if, on conclusion of test E517, the detectability criterion is not satisfied on any of the child blocks 
of the parent block B, the tested block B is held in memory in a memory-storage step E509 as a region available for 
55 inserting an information bit. 

[0396] In practice, a step of updating the capacity E51 0 can be implemented as described previously, in such a way 
that, for each block of coefficients tested which is suitable for insertion of an information bit, the value of the capacity 
is increased by one unit. 
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[0397] Hence, on conclusion of these adaptive partitioning steps E511 to E51 9, an image insertion capacity can be 

?™T' n ! ^ ,0 • 33 ° an 30 a PP r °P ria,e P*™™ °' the image corresponding to this capacity Q 

[0398] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit into the 

various regions of the partition. y 

[0399] Obviously, as before, the step E520 of insertion proper of the watermark can be coupled with the determination 

of the partition, ,n such a way that each block B available on conclusion of the detectability test E507, can actually be 

watermarked by the insertion of an information bit by modulation of the coefficients of this block B 

!^n?l L he ,nser,,on step E520 is carried out in 1he same way as when partition is achieved by iterative combination 

[0401] Moreover, in one enhanced embodiment of this technique of partition by iterative division, a division of a parent 

block B into child blocks can be adopted on completion of the detectability test E507 only if it represents a qain in 

capacitythat is to say into a quadtree-mode division only if the detectability criterion is satisfied over at least two child 

blocks. Otherwise, it may be preferable to keep the parent block of larger size so as to guarantee better invisibility of 

the inserted watermark ' 

[0402] Obviously, the partitioning techniques described previously are not in any way limiting 

[0403] Notably, it is possible to use various iterative segmentation techniques by division/merging, to which there 

belong the quadtree block partitioning described previously, but also graph partitioning. In general terms, there are the 

so-called bottom-up methods and the top-down methods depending on whether the iterative partitioning means is 

merging or division. " 

[0404] Hence it is possible, if appropriate, to combine two techniques of partitioning by division and combination 
This is the case when an initial partition is chosen, lor example, intermediate between the two extremes which are on 

1 !kT 3 Sln9 ' e re9i ° n andi ° P the ° ther hand ' a partition into re 9 ions °< the smallest admissible size, allowing 
both the operations of division and of combination. 

[0405] It is also possible to implement the method of determining a partition in order to maximise the capacity of the 
digital image by using several different segmentation techniques and then choosing the partition which confers the 
greatest capacity Q for the digital image. 

[0406] It is also possible to determine the best partition possible, that is to say the one yielding a set of blocks of 
coefficients of largest possible size, for inserting a fixed-length message, less than the maximum capacity of the digital 
image to be watermarked. a 
[0407] In all cases, the present invention makes it possible to generate and to insert a pseudo-random sequence 
once only over a set of coefficients which are available for carrying out the insertion, and, as a corollary, to apply only 
one single pre-distortion to the set of coefficients. 

[0408] It is thus possible rapidly to obtain a partition of the set of coefficients which makes it possible to insert a 
watermark with a fixed correct detection probability for an entire range of distortion lying between "no distortion" and 
predetermined distortion" applied during the determination of the partition. Moreover, the simulation of maximum dis- 
tortion before insertion is a quasi-deterministic operation which confers an advantage with respect to a theoretical 
statistical calculation. 

!!£?n] A description will be given now of the insertion method according to the invention in a preferred embodiment 
[O410] The watermarking signal is inserted after a transformation step E600, in a set of coefficients of the domain 
transformed into sub-bands. 

[0411] As illustrated in Figure 13, in a choosing step E601 , there is for example chosen the high-frequency sub-band 
r^V^r -rl resolution level ' corresponding to a high-pass filtering in the horizontal and vertical directions 
[041 2] There is thus a set of spatio-frequency coefficients, with a cardinal number equal to N, denoted for example 
X =( Xj, 1 < i < N}. r 

[0413] For an image I of size 512X512, the sub-band HH, is of size N = 256X256. 

[0414] Since the watermarking has to be imperceptible and indelible, and therefore difficult to locate and remove by 
pirating, ,t is chosen to .nsert a pseudo-random watermarking signal, spreading its spectrum in order to make this siqnal 
undetectable by spectral or statistical analysis. 

[0415] For example, a pseudo-random sequence w is considered, which follows a uniform law on the interval f-1 11 
wrth w =i Wj, 1 < i < PI, the length P of the sequence being less than or equal to N. 

[0416] Naturally, any pseudo-random watermarking signal, with a known distribution and a nil mean, may be suitable 
The most usual distributions for the watermarking signal w are, apart from the uniform distribution on [-1.1] mentioned 
above, the binary distribution i-1,1 , and the centred standardised Gaussian distribution N(0 1) 
[0417] The modulation which it is wished to apply to the coefficients in order to insert the watermarking signal uses 
a linear model of the type: 

X i ~ X i(i) + aw i' witn 1 - ) - p where 
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[0418] Xj (i) is a subset of spectral coefficients chosen from the set of coefficients X, and a is a weighting coefficient. 
[0419] Here, the weighting coefficient a is constant for any j. 

[0420] In this example embodiment, it is sought to determine whether the insertion of the watermarking signal w is 
actually possible in the set of coefficients X whilst meeting predetermined criteria of imperceptibility and detectability. 
5 [0421] According to the invention, the insertion method first of all includes a step E602 of calculating a maximum 
value of the weighting coefficient a as a function of the length P of the watermarking signal w ensuring imperceptibility 
of the watermarking signal w. 

[0422] This maximum value or visual amplitude o^, corresponds to the maximum modulation amplitude which can 
be used in the aforementioned linear insertion model beyond which an observer is capable of visually detecting a 

10 change at the level of the reconstructed watermarked image. 

[0423] In this embodiment, a visibility model is used which makes it possible to predict the visibility of a watermarking 
operation as a function of the different parameters, which are the representation of the signal through the spatio- 
I requency transformation T used, the sub-band Hh^ considered for the insertion, the type of distribution of the sequence 
of the watermarking signal w and the length P of the sequence w 

15 [0424] A simple model, developed by WATSON and described in the article entitled "Visibility of wavelet quantization 
noise", A.B. WATSON et al, IEEE Trans, on Image Process, 6(8), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of coefficients modulated from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for a detailed description of this model. 

[0425] During the step E602 of calculating the maximum value o^, use is thus made of a function dependent on the 
20 length P of the watermarking signal w and on the type of spectral transformation T used, but independent of the set of 
coefficients to be watermarked X. 
[0426] This calculation function used can be written: 

25 * V (P, S, w) = V Q base(S ! 1/p 

P 1/p (E[livr]) 

where 



30 <*base( s ) is a Dase value, dependent on the transformation used S and the (base) sub-band considered for the 

insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of this sub-band 
is modulated, 

P is strictly greater than 2, and 

£[lwl]P is the mathematical expectation of the function ItvlP 

35 

[0427] The base values a base (S) can be measured once and for all, for each sub-band of coefficients in the wavelet 
decomposition, from a single psychovisual measurement, and be stored in a table of visual amplitudes. 
[0428] p is the exponent of a Minkowsky sum and can be chosen equal to 5, for example. 
[0429] The mathematical expectation E[lwlP] corresponds to an estimation of the mean of the function itvlP. 
40 [0430] This visibility model does not take into account the image I to be watermarked itself and is independent thereof. 
This is equivalent to considering that the image I is uniform. It is a question of a "worst case" model since the presence 
of the image signal makes it possible to visually mask the modulation itself. 

[0431] In this example, this calculation step E602 is implemented in order to calculate the maximum value a v oi the 
weighting coefficient a tor a length P of the watermarking signal w equal to the cardinal number N of the set of modu- 
45 9 latable coefficients X. Thus a maximum value denoted cc^N) is calculated, which ensures imperceptibility of the wa- 
termarking signal w at the limit of perceptibility. 

[0432] In the example embodiment described here, for a value of p = 5, the visual amplitude for the set of coefficients 
of the sub-band Hh^ is given by: 

a v( N ) = ^02 



[0433] Next, in a calculation step E603, the length P of the watermarking signal w is calculated for a minimum value 
55 oq of the weighting coefficient a equal to said maximum value a^N) calculated and for a predetermined detection 
probability level. 

[0434] This is because, in order to meet a predetermined detectability criterion, the weighting coefficient a must be 
greater than a minimum value cxq referred to as the detection amplitude. 
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[0435] At the decoder., the received image I*, which corresponds to the watermarked initial image, possibly made 
noisy because of certain distortions applied to the image, is first of all transformed by a transformation identical to that 
used during the insertion of the watermarking signal in order to find once again the set of coefficients which are mod- 
ulated. 

[0436] ThusasetX^Xj*, 1 <i<N> is obtained, which represents the set of coefficients liable to have been modulated. 
By hypothesis, it is assumed that any distortion suffered by these coefficients can be modelled by an additive noise, 
decorrelated from the image signal itself n =*in jl 1 < i < NK. 
[0437] Each coefficient can thus be written: 

X-* = X- -»- aw- + n-, with 1 < i < N. 



[0438] The detection of the inserted watermarking signal w consists of constructing a hypothesis test for replying to 
the question "has the signal w been inserted in X*?". 
15 [0439] The advantage of constructing a hypothesis test is that its result makes it possible to give a theoretical con- 
fidence level to the decision since a hypothesis test is standardised with respect to known distributions. 
[0440] For a class of statistical tests based on correlation, it can be shown that this detection amplitude can be 
calculated from a calculation function of the type: 



aQo 2 x +o 2 n ) 



where a* is equal to the variance of the coefficients to be modulated, is equal to the variance of the additive 
25 noise and a, b and c are constants which depend on the distribution of the watermarking signal w and on the required 
detection probability level p. 

[0441] Thus this minimum value ocq is determined according notably to the length P of the watermarking signal w, a 
detection probability level p, the variance of the coefficients to be modulated, the distribution of the watermarking signal 
and possibly the variance of an additive noise decorrelated from the coefficients modelling any distortions applied to 
30 the coefficients. 

[0442] Conversely, it is also possible to express the minimum detection length, denoted P°, as a function of the 
weighting coefficient a, for a fixed detection probability level: 

2 2 2 2 

35 P°(a) = — f 

bo. 

[0443] It would be possible to use this result to check the parameters which make it possible to insert a watermarking 
signal w which will be robust to any additive noise with an energy less than or equal to a n 2 . 

40 [0444] As illustrated on the curves in Figure 14, from the visual amplitude cty calculated as a function of the length 
P of the watermarking signal to be inserted, the insertion method includes, in general terms, a step of determining 
compatible values of the length P of the watermarking signal w, of the weighting coefficient a and of the probability 
level p so that the value of the weighting coefficient determined a is less than or equal to the maximum value or visual 
amplitude a v and greater than or equal to the minimum value or detection amplitude a D for these determined values 

45 of the length P of the watermarking signal w and the detection probability level p. 

[0445] The weighting coefficient a, the length P of the watermarking signal w and the detection probability level p 
must thus be determined so that the following inequality is satisfied: 

so a D < a < a v 

[0446] As illustrated in Figure 14, in order to calculate the minimum value P° of the watermarking signal w for which 
the aforementioned inequality will be satisfied, the formula linking the length P of the watermarking signal w to the 
amplitude of the weighting coefficient a is used, with a = a^N) and w ' tn a detection probability level of p%. 
55 [0447] In this example a standardised statistical test is taken, based on the correlation between the watermarked 
signal X* and the inserted watermarking signal w, as proposed for example in the document entitled "On Resolving 
Rightful Ownerships of Digital Images by Invisible Watermarks" by W. ZENG and B. LIU, in Proc. ICIP 97, pages 
552-555, October 1997. 
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[0448] It can be shown that, for this test: 

[0449] A calculation step E604 thus makes it possible to calculate the three aforementioned constants a, b and c so 
that 

a= t 0 JE(w 2 ) t 

and 




[0450] The values £( vfi) and o 2 2 can be calculated theoretically as a function of the distribution of w. For the uniform 
distribution considered here, theTheoretical values are E(vfi)= 1/3 and o 2 = 4/45 

[0451] The value ^ is a threshold value, dependent on the probability oT correct detection p:t 0 = 2xi p . The values 
of t p are percentile of a centred Gaussian law and have a variance equal to 1 . These values are tabulated in the state 
of the art, for example in the book entitled "Probability and Statistics" by A. PAPOULIS, Prentice-Hall, 1990. 
[0452] For example, t p = 3,291 corresponds to p = 0,9995, that is to say to a detection probability level equal to 

99,95%. 

[0453] In this numerical example it is considered that any additive noise is not taken into account so that o n = 6. 
[0454] In an estimation step E605, the variance of the set of coefficients X is also estimated. The variance is estimated 
according to the following formula: 



<x 3 v=— f{X ( -M x Y 

with 



1 N 

[0455] Then, for a detection probability fixed at 99.95%, the function linking the length P of the watermarking signal 
w to be inserted and the weighting coefficient a is then derived therefrom: 

14,44 X 0%+ 3,34a 2 

P= S 

0,11a 

[0456] At the calculation step E603, the minimum length P° is determined from the above formula for a weighting 
coefficient equal to the maximum value a v (N) calculated at the calculation step E602. 

[0457] The insertion method includes, in this example embodiment, a step E606 of comparing this calculated length 
P° with the cardinal number N of the set of modulatable coefficients X. 

[0458] If, at this comparison step E606, the calculated length P° is less than or equal to the cardinal number N of 
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the set of modulatable coefficients X, it is deduced therefrom that the insertion of a watermarking signal w of length P° 
is possible and meets the required imperceptibility and detection criteria. 

[0459] In reality, it is possible to effect an insertion of a watermarking signal w which meets these criteria since the 
length P of this watermarking signal lies between the calculated value P° and the cardinal number N of the set of 
modulatable coefficients X. 

[0460] Then the insertion proper of a watermarking signal of length P is carried out by modulating a subset of coef- 
ficients with a cardinal number P according to the linear modulation level, with a weighting coefficient a equal to the 
calculated maximum value <x v (N). 

[0461] The insertion method includes, in the usual fashion, a reverse spatio-frequency transformation step which 
makes it possible to obtain, at the output of the coder, the watermarked image I*. 

[0462] Where insertion is possible, it is also possible to calculate the number of theoretical bits which can be inserted. 
This is because, at the decoder, the decoding operation consists of detecting whether a pseudo-random sequence w 
has been inserted. The response being binary (yes/no), it can be considered that a pseudo-random sequence of length 
P makes it possible to insert and extract an information bit in the image. By repeating the insertion operation on several 
subsets of spectral coefficients, it is thus possible to insert and extract several information bits. 
[0463] Consequently, since in this example there are a set of coefficients of cardinal number N and since only a 
number P° of coefficients are necessary to insert and detect an information bit, it is possible to derive therefrom the 
total number Q of information bits which can be inserted, referred to as the capacity of the image: Q = N/P°. 
[0464] Moreover, if at the comparison step E606 the calculated length P° is greater than the cardinal number N of 
the set of modulatable coefficients X, it is a priori deduced therefrom that the insertion of a watermarking signal meeting 
the criteria of imperceptibility and detectability is not possible. 

[0465] In this case, the insertion method can if necessary include an additional step of calculating the detection 
probability level p' obtained for a minimum value a D of the weighting coefficient equal to the maximum value ct v (N) 
calculated for a length P of the watermarking signal w equal to the cardinal number N of the set of modulatable coef- 
ficients X. 

[0466] If this detection probability level, which will be less than that fixed initially at 99.95%, is judged satisfactory, 
the insertion of a watermarking signal of length equal to N can be effected, that is to say all the coefficients of the sub- 
band of frequency HH^ are modulated. 

[0467] In a second embodiment of the invention, the insertion method can include a step of calculating a threshold 
value P* of the length of the watermarking signal w, as illustrated in Figure 14. 

[0468] This threshold value P* is determined so that the minimum value a D of the weighting coefficient is equal to 
the maximum value cc v of the weighting coefficient for a length P of the watermarking signal w equal to this threshold 
value P* and a predetermined detection probability level p. 

[0469] This calculation step makes it possible to determine, from formulae expressing the detection amplitude a D 
and the visual amplitude a v notably as a function of the length P of the watermarking signal w, the minimum value P* 
of the insertion signal from which the weighting coefficient a can satisfy the inequality a D < a <a v . 
[0470] In Figure 14, this minimum value P* is given by the intersection of the two curves representing the variation 
in the visual amplitude a v and the variation in the detection amplitude as a function of the length P of the inserted 
watermarking signal w, for a predetermined detection probability level. 

[0471] Thus, in general terms, the insertion method according to the invention has the advantage of being able to 
control a priori and theoretically the detection probability and imperceptibility of the watermarking signal which it is 
sought to insert in an image. 

[0472] Naturally, many modifications could be made to the example embodiments described above without departing 
from the scope of the invention. 

[0473] Thus the insertion technique used could be applied to the raw digital image, without undergoing any spatio- 
frequency transformation prior to the modulation of the coefficients. 

[0474] In this case, the modulated coefficients are coefficients representing the digital image solely in the spatial 
domain. 

[0475] Moreover, the spatio-frequency transformation applied to the image can use analysis and synthesis filters 
other than those described above, or even be a transformation other than discrete Fourier transformation, discrete 
cosine transformation or Fourier-Mellin transformation. These transformations are used currently in conventional 
processings of digital images. 

[0476] In addition, the invention is not limited to the example embodiments described but also concerns any embod- 
iment where the values of the weighting coefficient a, of the length P of the watermarking signal and of the detection 
probability level p are determined in a correlated fashion so as to satisfy the inequality Oq < a <ay. 
[0477] The present invention also concerns a device for determining a partitioning and a decoding device which are 
adapted to implement the partitioning and decoding methods described previously. 

[0478] These devices for determining a partitioning and decoding can be used in a computer 50 as illustrated in 
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Figure 15, independently or in the same computer 50. 

[0479] The partitioning determination device comprises partitioning means 500, 502, 503 adapted to perform an 
adaptive partitioning according to a criterion of detectability of an information bit inserted on each region. 
[0480] These partitioning means use on coding a maximum value a v of the weighting coefficient ensuring impercep- 
s tibility of the watermarking signal for modulating all the coefficients representing said digital image, possibly with a 
modulation factor kj dependent on each coefficient, as described previously. 

[0481] This device for determining a partitioning also has means 500, 502, 503 of comparing the size P of each 
region of the partitioning with a minimum size P s1al corresponding to minimum size of a statistically significant sample 
for the detection of a watermarking signal inserted in said region. 
10 [0482] The device also has means 500, 502, 503 for applying a pre -distort ion before calculating a detectability cri- 
terion. 

[0483] The partitioning means comprise, means 500, 502, 503 of modulating all the coefficients of each region by 
inserting the watermarking signal, means 500, 502, 503 of calculating a detectability amplitude from the detectability 
criterion and means 500, 502, 503 of validating the detectability criterion, adapted tocompare the detectability amplitude 

15 with a predetermined threshold value Tc. 

[0484] The partitioning means 500, 502, 503 are adapted to effect a partitioning by iterative division of the numerical 
image and, for each region at a rank in the said partitioning, the validation means 500, 502, 503 are adapted to validate 
a partitioning of a rank higher in the region if and only if there exist at least two subregions of the region for which the 
detectability criterion is validated. 

20 [0485] This device also comprises means 500, 502, 503 of incrementing the value of a capacity Q of the digital image, 
adapted, for each region at a rank in the partitioning, to increment the capacity value Q when the partitioning of rank 
directly higher is not validated, the capacity Q being equal to the cardinal number of the distinct regions obtained by 
the partitioning on which the detectability criterion- is validated. 

[0486] In another embodiment: this device for determining a partitioning comprises incrementation means 500, 502, 
25 503 adapted, for each region R at a rank in the partitioning, to increment the value of a capacity Q when the partitioning 
of rank directly higher is not validated, comparison means 500, 502, 503 adapted to compare the value of the capacity 
Q with a fixed capacity value Qs and reiteration means 500, 502, 503 adapted to reiterate the partitioning for a region 
having the highest detectability amongst all the other regions to be processed when the capacity value Q is less than 
the fixed capacity value Qs. 

30 [0487] In a third embodiment of the invention, the partitioning means 500, 502 : 503 are adapted to effect a partitioning 
by iterative merging of the digital image and there are adapted, for each region of the partitioning for which the de- 
tectability criterion is not validated, to merge the region with another region of the partitioning, and preferably with 
another region of the partitioning for which the detectability criterion is not validated. 

[0488] In an entirely similar fashion, the device for decoding a watermarking signal in a digital image comprises 
35 means 500, 502, 503 of partitioning the digital image to be decoded into distinct regions, adapted to perform an adaptive 
partitioning as a function of a criterion of detectability of an information bit obtained by demodulating the coefficients 
on each region. 

[0489] It also has means 500, 502, 503 of comparing the size P of each region R of the partitioning with a minimum 
size P stat corresponding to the minimum size of a statistically significant sample for the detection of a watermarking 
40 signal inserted in the region. 

[0490] It comprises means 500, 502, 503 of calculating a detection amplitude from the detectability criterion and 
means 500, 502, 503 of validating the detectability criterion, adapted to compare the detection amplitude with a pre- 
determined threshold value Td for decoding. 

[0491] This predetermined threshold value Td for decoding is less than the predetermined threshold value Tc used 
45 during the partitioning determination method described previously. 

[0492] In a preferred embodiment, the device for determining a partition includes means 500, 502, 503 for checking 
a criterion of detectability, over each region, of an information bit inserted by means of a pseudo-random sequence 
with average of zero by modulation of the coefficients of the region. " ■ 

[0493] It also includes means 500, 502, 503 for generating a centred pseudo-random sequence equal in size to the 
50 cardinal number of the set of coefficients representative of an image, here equal to N. 

[0494] Means 500, 502, 503 of modulating the set of coefficients X by the centred pseudo-random sequence in order 
to insert the same information bit onto this set of coefficients are also incorporated into the device for determining a 
partition 11 as shown on figure 1. 

[0495] The latter may also includes means 500 : 502, 503 for applying a distortion D which makes it possible to 
55 provide for and guarantee the robustness of the inserted watermark in the case of certain types of distortion, typically 
a compression of the digital image in order to store it. 

[0496] These means of applying a distortion make it possible actually to simulate a predetermined distortion over 
the set of modulated coefficients. Here they comprise, whenever the coefficients of the set X are spectral coefficients, 
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means 500, 502, 503 for inverse spectral transformation of the set of coefficients modulated in order to obtain a wa- 
termarked image : means 500, 502, 503 for applying a distortion to this watermarked image, and means 500, 502, 503 
for spectral transformation of the watermarked image in order to recover the set of modulated coefficients after distor- 
tion. 

[0497] In this example, the inverse-spectral-transformation means are suitable for implementing multi-resolution 
spectral wavelet recomposition and the forward spectral transformation means are suitable for again reiterating the 
breakdown of the digital image I into wavelets. 

[0498] The device for determining a partition 11 , in a first way. also includes means for partitioning by iterative com- 
bination, on the basis of an initial segmentation of the set of coefficients X. 

[0499] For preference, this initial segmentation makes it possible to segment the set of coefficients X into a set of 
regions of a size which is statistically significant for validating the detectability criterion over these regions. 
[0500] The means of partitioning 500, 502, 503 by iterative combination are suitable, for each region of the partition 
for which the detectability criterion is not satisfied, for combining the region with another region of the partition. 
[0501] Obviously, the partitioning means 500, 502, 503 could also be suitable for forming this partition of the image 
by iterative division of the digital image according to a second way of performing the preferred embodiment of the 
invention. 

[0502] This partition-determining device 11 makes it possible, as output, to obtain a partition into regions of the set 
of coefficients X, with it being possible to insert an information bit onto each of these regions with a reliable detection 
probability and robustness to a given distortion. The maximum number of regions which can be obtained by this adaptive 
partitioning of the set X corresponds to the practical insertion capacity of this set X. 

[0503] As illustrated in Figure 1 , an insertion device 10, incorporating the device for determining a partition 11 de- 
scribed above, moreover comprises means 13 for reading a watermark to be inserted, containing several information 
bits. By way of an example which is not at all limiting, a binary watermark may be 1000111 , representing, for example, 
the coding of the name of the author of the image to be watermarked. 

[0504] The insertion means 1 2 proper may be conventional and make it possible, by modulation, to insert the various 
information bits of the watermark onto the various insertion carriers determined by the device for determining a partition 



[0505] According to the invention, the device for inserting a watermarking signal w comprises means 500, 502, 503 
of spatio-frequency transformation of an image I, and for example analysis filters associated with decimators by two 
adapted to effect a wavelet decomposition of an image I It also has reverse spatio-frequency recomposition means 
500, 502, 503 for recomposing the image I after the insertion of the watermarking signal in the domain transformed 
into sub-bands. 1 

[0506] It also comprises means 500, 502, 503 of inserting a watermarking signal w of length P adapted to modulate 
a subset of coefficients of cardinal number P according to the linear model described previously using a weighting 
coefficient a. 

[0507] This insertion device also comprises means 500, 502, 503 of calculating the maximum value cc v of the weight- 
ing coefficient as a function of the length P of the watermarking signal w ensuring imperceptibility of the watermarking 
signal. These calculation means use the function a v (P, T) described above and dependent on the length P of the 
watermarking signal w and the spectral transformation T used for decomposing the image signal I. 
[0508] in general terms, this insertion device also has means 500, 502, 503 of determining compatible values of the 
length P of the watermarking signal w, the weighting coefficient a and the probability level p so that the value of the 
weighting coefficient determined is less than or equal to the maximum value and greater than or equal to the minimum 
value Oq for the determined values of the length P of the watermarking signal w and the detection probability level p. 
[0509] The calculation means 500, 502, 503 are also adapted to determine a minimum value or detection amplitude 
a D using the calculation functioning described above of the type: 



[0510] In one embodiment of the invention, the insertion device also has means of calculating 500 : 502, 503 the 
length P° of the watermarking signal wfor a minimum value a D of the weighting coefficient equal to a maximum value 
a v (N) calculated for a length P of the watermarking signals w equal to the cardinal number N of the set of modulatable 
coefficients ensuring imperceptibility of the watermarking signal w at the limit of perceptibility and for a predetermined 
detection probability level p; and means 500, 502, 503 of comparing this length P° calculated with the cardinal number 
N of the set of modulatable coefficients X. 

[0511] It also comprises means 500, 502, 503 of calculating a detection probability level p' obtained for a minimum 
value Oo of the weighting coefficient equal to the maximum value a v calculated for a length P of the watermarking 



11. 



« D (P) = 
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signal equal to the cardinal number N of the set of modulatable coefficients X. 

[0512] In another embodiment, the insertion device comprises means 500, 502, 503 of calculating a threshold value 
P* of the length of the watermarking signal w determined so that the minimum value cio of the weighting coefficient is 
equal to the maximum value a v of the weighting coefficient of the length P of the watermarking signal w equal to this 

5 threshold value P* and a predetermined detection probability level p. 

[0513] All the previously described devices are incorporated in a microprocessor 500 of the computer 50, a read- 
only memory 502 being adapted to store a program for determining a partitioning and/or inserting a watermarking 
signal and/or decoding the watermarking signal and a random access memory 503 comprising registers adapted to 
store variables modified during the running of the programs. 

10 [0514] The microprocessor 500 is integrated into the computer 50, which can be connected to different peripherals, 
for example a digital camera 507 or a microphone 511, by means of an input/output card 506 in order to receive and 
store documents. 

[0515] The digital camera 507 makes it possible notably to supply images to be authenticated by inserting a water- 
marking signal. 

15 [0516] This computer 50 has a communication interface 512 connected to a communication network 51 3 in order to 
receive, if necessary, images to be watermarked. 

[0517] The computer 50 also has document storage means, such as a hard disk 508, or is adapted to cooperate, by 
means of a disk drive 509, with removable document storage means such as diskettes 510. 

[0518] These fixed or removable storage means can also contain the code of the different methods according to the 
20 invention which, once read by the microprocessor 500, will be stored on the hard disk 508. 

[0519] By way of variant, the program enabling the devices to implement the invention can be stored in the read only 
memory 502. 

[0520] As a second variant, the program can be received in order to be stored as described above by means of the 
communication network 513. 

2S [0521] The computer 50 also has a screen 504 for serving, for example, as an interface with an operator by means 
of the keyboard 51 4 or any other means. 

[0522] The central unit 500 (CPU) will execute the instructions relating to the implementation of the invention. On 
powering up, the programs and methods relating to the invention stored with a non-volatile memory, for example the 
read only memory 502, are transferred into the random access memory 503 (RAM) which will then contain the exe- 

30 cutable code of the invention and the variables necessary for implementing the invention. 

[0523] This random access memory 503 contains different registers for storing the variables necessary for the running 
of the program, and notably a register for storing the coefficients of the regions R which are temporarily modulated or 
demodulated in order to verify the detectability criterion, the minimum value P stat , and the threshold value Tc and Td 
respectively on coding and decoding. 

3S [0524] This random-access memory 503 also includes an Xj, Xj register for storing the coefficients of the regions 
which are modulated or demodulated temporarily so as to check the detectability criterion and a register w for storing 
the centred pseudo-random sequence generated in order to determine the partition of the set X. 
[0525] Finally, the random access memory 503 contains a register for storing the table of visual amplitudes of a base 
(T), a register for storing the spectral coefficients X, a register for calculating the visual amplitude a v (N) at the limit of 

40 imperceptibility, a register for storing the values of the constants a, b, c and of the variance a x 2 necessary for determining 
the detection amplitude a D as a function of the length P of the watermarking signal w and a register for storing the 
minimum value P° calculated. 

[0526] A communication bus 501 affords communication between the different sub-elements of the computer 50 or 
linked to it. The representation of the bus 501 is not limitative and notably the microprocessor 500 is able to indicate 
45 instructions to any sub-element directly or by means of another sub-element. 



Claims 

50 1. Method of determining a partitioning into distinct regions of a digital image (I) in order to insert a watermarking 
signal (W) by modulating coefficients (X) of said regions (R) representing a digital image (I), characterised in that 
the partitioning into regions is effected by an adaptive partitioning as a function of a criterion of detectability of an 
information bit inserted on each region (R). 

55 2. Method of determining a partitioning according to Claim 1 , characterised in that the detectability criterion is calcu- 
lated by using a maximum value (a v ) of the weighting coefficient ensuring imperceptibility of the watermarking 
signal for modulating all the coefficients representing said digital image. 
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3. Method of determining a partitioning according to claim 2, characterised in that a maximum value (a v ) can be 
calculated by: 



a v (P,S, iv) = base 



P 1/p (E[l^ P ]) 1/,i 



where 

10 <*base(S) ' s a base value, dependent on the transformation used S and on the (base) sub-band considered for 

the insertion, of the maximum weighting coefficient ensuring imperceptibly when a single coefficient of said 
sub-band is modulated, 
P is strictly greater than 2, and 

£[lwlP] is the mathematical expectation of the function IwlP. 

15 

4. Method of determining a partitioning according to one of Claims 2 or 3 : characterised in that the detectability 
criterion is calculated by using, tor each coefficient (Xj) to be modulated, a weighting coefficient in accordance with 
a law of the form a t = kj.a v , where kj is a modulation factor dependent on the coefficients situated close to the 
relevant coefficient (Xj) on said region and is equal to said maximum value of the weighting coefficient. 

20 

5. Method of determining a partitioning according to one of claims 1 to 4, characterised in that it comprises a step 
(E307) of applying a distortion before calculating the detectability criterion. 

6. Method of determining a partitioning according to one of Claims 1 to 5, characterised in that it includes a prior step 
25 (E105) of comparing the size (P) of each region (R) of the partitioning with a minimum size (P stat ) corresponding 

to the minimum size with a statistically significant sample for the detection of a watermarking signal inserted in 
said region. 

7. Method of determining a partitioning according to one of Claims 1 to 6, characterised in that, for each region (R) 
30 of the partitioning, all the coefficients of said region are modulated by inserting the watermarking signal (W), a 

detectability amplitude (T) is calculated from said detectability criterion and said detectability criterion is validated 
' by comparing said detectability amplitude (T) with a predetermined threshold value (Tc). 

8. Method of determining a partitioning according to claim 7, characterised in that said detectability criterion is vali- 
35 dated if a minimum modulation length (P min (x)) calculated for ensuring detectability of the inserted watermarking 

signal (W) is less than or equal to the cardinal number of said region (R). 

9. Method of determining a partitioning according to claim 1 , characterised in that it comprises the following steps: 

40 - generation (E501 ) of a centred pseudo-random sequence (w) equal in size to the cardinal number (N) of a set 

of coefficients representative of at least a part of an image (I), formed from centred pseudo-random sub- 
sequences (w*); 

modulation (E502) of the said set of coefficients by the said centred pseudo-random sequence (w) in order to 
insert the same information bit on the said set of coefficients; and 
45 - checking (E506, E507, E514, E516), over each region, of a criterion of detectability of the information bit 

inserted by modulation. 

10. Method of determining a partitioning according to Claim 9, characterised in that it comprises a single step of applying 
a distortion (E503) to the said set of modulated coefficients, before the step of checking (E506, E507, E514, E516) 

so of a detectability criterion. 

11. Method of determining a partitioning according to Claim 10, characterised in that the said set of coefficients is a 
set of spectral coefficients (X) originating from a spectral transformation of a digital image, and in that the said step 
(E503) of applying a distortion includes the following sub-steps: 



inverse spectral transformation (E503a) of the said set of modulated coefficients in order to obtain a water- 
marked image; 

application (E503b) of a distortion to the said watermarked image; and 
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spectral transformation (E503c) of the watermarked image in order to obtain a set of modulated coefficients 
after distortion. 

12. Method of determining a partitioning according to one of Claims 9 to 11 , characterised in that the centred pseudo- 
5 random sequence (w) is formed from pseudo-random sub-sequences (wk) of identical size. 

13. Method of determining a partitioning according to one of Claims 9 to 1 2, characterised in that the pseudo-random 
sub-sequences (wk) have a size greater than or equal to a minimum number of coefficients forming a region of 
the said partition of a size which is statistically significant for validating the said criterion of detectability over the 

io said region. 

14. Method of determining a partitioning according to one of Claims 9 to 13, characterised in that, at the generation 
step (E501 ), each centred pseudo-random sub-sequence (wk) is created by generating half of the pseudo-random 
samples of the said sub-sequence via a pseudo-random number generator and by symmetrising the said pseudo- 

is random samples generated. 

15. Method of determining a partitioning according to Claim 14, characterised in that, at the generation step (E501), 
each centred pseudo-random sub-sequence (wk) is created by moreover permutating the said pseudo-random 
samples using a secret key. 

20 

16. Method of determining a partitioning according to one of Claims 9 to 13, characterised in that, at the generation 
step (E501), each centred pseudo-random sub-sequence (wk) is created by generating pseudo-random random 
samples via a pseudo-random number generator and by redistributing over each sample the sum of the said 
pseudo-random samples generated. 

25 

17. Method of determining a partitioning according to one of Claims 7 to 16, characterised in that a partitioning by 
iterative division (E103) of the digital image is effected and in that, for each region at a rank of said partitioning, a 
partitioning of a directly higher rank of said region is validated if and only if there exist at least two sub-regions of 
said region for which said detectability criterion is validated. 

30 

18. Method of determining a partitioning according to Claim 17, characterised in that, for each region at a rank of said 
partitioning, the value of a capacity (Q) of the digital image is incremented (E112) when the partitioning of directly 
higher rank is not validated, said capacity (Q) being equal to the cardinal number of the distinct regions obtained 
by said partitioning on which the detectability criterion is validated. 

35 

19. Method of determining a partitioning according to Claim 18, characterised in that, for each region at a rank of said 
partitioning, the value of a capacity (Q) is incremented (El 12) when the partitioning of immediately high rank is 
not validated, said value of the capacity (Q) is compared (E11 5) with a fixed capacity value (Qs) and the partitioning 
is reiterated for a region having the highest detectability amplitude amongst all the other regions to be processed 

40 when said capacity value is lower than said fixed capacity value. 

20. Method of determining a partitioning according to one of Claims 7 to 16, characterised in that a partitioning by 
iterative merging of the digital image is effected and in that, for each region of said partitioning for which said 
detectability criterion is not validated, said region is merged with another region of said partitioning. 

45 

21 . Method of determining a partitioning according to Claim 20, characterised in that said other region of the partitioning 
is if possible a region for which said detectability criterion is not validated. 

22. Method of inserting, into a digital image, a watermark comprising one or more information bits, characterised in 
so that it comprises a method of determining a partition according to one of Claims 1 to 21, and a step (E520) of 

inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 

23. Method of decoding a watermarking signal in a digital image, inserted by modulating coefficients (X) representing 
said image in distinct regions forming a partitioning of said digital image, characterised in that a partitioning into 

55 distinct regions of the digital image to be decoded is determined by effecting an adaptive partitioning as a function 

of a criterion of detectability of an information bit obtained by demodulating coefficients on each region (R). 

24. Decoding method according to Claim 23, characterised in that it includes a prior step (E201 ) of comparing the size 
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(P) of each region of the partitioning with a minimum size (P slal ) corresponding to the minimum size of a statistically 
significant sample for the detection of a watermarking signal inserted in said region. 

25. Decoding method according to one of Claims 23 or 24, characterised in that, for each region (R) of the partitioning. 
s a detection amplitude (T) is calculated (E203) from said detectability criterion and said detectability criterion is 

validated (E203) by comparing said detection amplitude (T) with a predetermined threshold value (Td) for decoding. 

26. Decoding method according to Claim 25, characterised in that the predetermined threshold value (Td) for decoding 
is less than said predetermined threshold value (Tc) used during the partitioning determination method according 

io to Claim 7. 

27. Device for determining a partitioning into distinct regions of a digital image for inserting a watermarking signal by 
modulating coefficients of said regions representing a digital image, characterised in that it comprises partitioning 
means (500, 502, 503) adapted to effect an adaptive partitioning as a function of a criterion of detectability of an 

*5 information bit inserted on each region. 

28. Partitioning determination device according to claim 27, characterised in that it also includes means (500, 502, 
503) for applying a distortion before calculating the detectability criterion. 

20 29. Partitioning determination device according to one of claims 27 to 28, characterised in that it also has means (500, 
502, 503) of comparing the size of each region of the partitioning with a minimum size corresponding to the minimum 
size of a statistically significant sample for the detection of a watermarking signal inserted in the said region. 

30. Partitioning determination device according to one of Claims 27 to 29, characterised in that the partitioning. means 
25 (500, 502, 503) comprise means (500, 502, 503) of modulating all the coefficients of each region by inserting the 

watermarking signal, means (500, 502, 503) of calculating a detectability amplitude from said detectability criterion 
and means (500, 502, 503) of validating said detectability criterion, adapted to compare said detectability amplitude 
with a predetermined threshold value 

30 31. Partitioning determination device according to claim 27, characterised in that it includes: 

means (500, 502, 503) for generating a centred pseudo-random sequence (w) equal in size to the cardinal 
number (N) of a set of coefficients representative of at least a part of an image (I), formed by centred pseudo- 
random sub-sequences (w k ); 

35 - means (500, 502, 503) for modulation of the said set of coefficients by the said centred pseudo-random se- 

quence (w) in order to insert the same information bit on the said set of coefficients; and 
means 500, 502, 503) for checking, over each region, a criterion of detectability of the information bit inserted 
by modulation. 

40 32. Partitioning determination device according to Claim 31 , characterised in that it comprises means (500, 502, 503) 
for applying a distortion to the said set of modulated coefficients. 

33. Partitioning determination device according to Claim 32, characterised in that the said set of coefficients is a set 
of spectral coefficients (Xj) originating from a spectral transformation of a digital image, and in that the said means 

45 for applying a distortion comprise: 

means (500, 502, 503) for inverse spectral transformation of the said set of modulated coefficients in order to 
obtain a watermarked image; 

means (500, 502, 503) for application of a distortion to the said watermarked image; and 
50 - means (500, 502, 503) for spectral transformation of the watermarked image in order to obtain a set of mod- 

ulated coefficients after distortion. 

34. Partitioning determination device according to one of Claims 30 to 33, characterised in that the partitioning means 
(500, 502 : 503) are adapted to effect a partitioning by iterative division of the digital image and in that, for each 

55 region at a rank of said partitioning, the validation means are adapted to validate a partitioning of a directly higher 

rank of said region if and only if there exist at least two sub-regions of said regions for which said detectability 
criterion is validated. 
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35. Partitioning determination device according to Claim 34, characterised in that it also comprises means (500, 502. 
503) of incrementing the value of a capacity of the digital image adapted, for each region at a rank of said parti- 
tioning, to increment the said capacity value when the partitioning of directly higher rank is not validated, said 
capacity being equal to the cardinal number of the distinct regions obtained by said partitioning on which the 

5 detectability criterion is validated. 

36. Partitioning determination device according to Claim 34, characterised in that it also comprises incrementation 
means (500, 502, 503) adapted, for each region at a rank of said partitioning, to increment the value of a capacity 
when the partitioning of directly higher rank is not validated, comparison means (500, 502, 503) adapted to compare 

io said value of the capacity with a fixed capacity value and reiteration means (500, 502, 503) adapted to reiterate 

the partitioning for a region having the highest detectability amplitude amongst all the other regions to be processed 
when said capacity value is lower than said fixed capacity value. 

37. Partitioning determination device according to one of Claims 30 to 33, characterised in that the partitioning means 
is (500, 502, 503) are adapted to effect a partitioning by iterative merging of the digital image and are adapted, for 

each region of said partitioning for which said detectability criterion is not validated, to merge said region with 
another region of said partitioning. 

38. Partitioning determination device according to Claim 37, characterised in that said other region of the partitioning 
20 is if possible a region for which said detectability criterion is not validated. 

39. Partitioning determination device according to one of Claims 27 to 38, characterised in that it is incorporated in a 
microprocessor (500), a read only memory (502) being adapted to store a program for determining a partitioning 
and a random access memory (503) comprising registers adapted to store variables modified during the running 

2S of said program. 

40. Device for inserting, into a digital image, a watermark comprising one or more information bits, characterised in 
that it comprises a device (1 1 ) for determining a partition according to one of Claims 27 to 39, and means (12) for 
inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 

30 

41 . Insertion device according to Claim 40, characterised in that it is incorporated into a microprocessor (500), a read- 
only memory (502) being suitable for storing in memory a program for inserting a watermark by modulation, and 
a random-access memory (503) including registers suitable for storing in memory the variables modified during 
the running of the said program. 

35 

42. Device for decoding a watermarking signal in a digital image, inserted by modulation of coefficients representing 
said image in distinct regions forming a partitioning of said digital image, characterised in that it comprises means 
(500, 502, 503) of partitioning of the digital image to be decoded into distinct regions, adapted to effect an adaptive 
partitioning as a function of a criterion of detectability of an information bit obtained by demodulation of the coef- 

40 ficients on each region. 

43. Decoding device according to Claim 42, characterised in that it also has means (500, 502, 503) of comparing the 
size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistical 
significant sample tor the detection of a watermarking signal inserted in said region. 

45 

44. Decoding device according to one of Claims 42 or 43, characterised in that it comprises means (500, 502, 503) of 
calculating a detection amplitude from said detectability criterion and means (500, 502, 503) of validating said 
detectability criterion adapted to compare said detection amplitude with a threshold value predetermined for the 
decoding. 

so 

45. Decoding device according to Claim 44, characterised in that the predetermined threshold value for the decoding 
is less than said predetermined threshold value used during the partitioning determination method according to 
Claim 7. 

s & 46. Decoding device according to one of Claims 42 to 45, characterised in that it is incorporated in a microprocessor 
(500), a read only memory (502) being adapted to store a program for decoding a watermarking signal and a 
random access memory (503) comprising registers adapted to store variables modified during the running of said 
program. 
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47. Method of processing an image, characterised in that it comprises the following steps: 

inputting an image (E500); 

inserting a first watermarking information into said inputted image (E502); 
applying a distortion to the image obtained in said insertion step (E503a-c); 

detecting said first watermarking information from the image obtained in said distortion applying step 
(E504-E508); and 

determining an amount of a second watermarking information which can be inserted into said inputted image 
in accordance with the result of the detection obtained in said detection step (E510). 

48. Processing method according to Claim 47, characterised in that it futher comprises the step of displaying the 
amount of the watermarking information determined in said determining step. 

49. Image processing apparatus, characterised in that it comprises: 



means (500, 502, 503) for inputting an image; 

means (500, 502, 503) for inserting a first watermarking information into said inputted image; 
means (500, 502, 503) for applying a distortion to the image obtained by said insertion means; 
means (500, 502, 503) for detecting said first watermarking information from the image obtained by said dis- 
20 lortion applying means; and 

means (500, 502, 503) for determining an amount of a second watermarking information which can be inserted 
into said inputted image in accordance with the result of the detection obtained by said detection means. 



50. Image processing apparatus according to Claim 49, characterised in that it further comprises the means (500, 502, 
503) for displaying the amount of the watermarking information determined by said determining means. 



51. Method of inserting a watermarking signal (w) in a set of coefficients (X) representing a digital image (I), in which 
at least one subset of coefficients is modulated according to a linear model using a weighting coefficient (a), the 
weighting coefficient (a) being greater than a minimum value (a D ) determined as a function notably of the length 
30 (P) of the watermarking signal (w) and a detection probability level (p), characterised in that it includes: 

a step (E602) of calculating a maximum value (a v (N)) of the weighting coefficients functional to the length (P) 
of the watermarking signal (w) ensuring imperceptibility of the watermarking signal. 

3S 52. Watermarking signal insertion method according to Claim 51 , characterised in that it also includes: 

a step (E603) of determining compatible values of the length (P) of the watermarking signal (w) of the weighting 
coefficient (a) and of the detection probability level (p) so that the value of the weighting coefficient determined 
(a) is less than or equal to said maximum value (a v ) and greater than or equal to said minimum value (<x D ) for 
said given values of the length (P) of the watermarking signal (w) and the detection probability level (p). 

53. Watermarking signal insertion method according to one of Claims 51 or 52, characterised in that, during the step 
(E602) of calculating the said maximum value (a v ), a function is used which is independent of said set of coefficients 
to be watermarked (X). 

54. Watermarking signal insertion method according to one of Claims 51 to 53, characterised in that said coefficients 
are spatio-frequency coefficients obtained by spatio-frequency transformation of an image (I), the calculation of 
said maximum value (ocj of the weighting coefficient depending on the type of transformation used (S). 

so 55. Watermarking signal insertion method according to Claim 54, characterised in that said calculation function used 
is of the type 



v ' ' ' p ,/p (£[iwiy p 

where 
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a base (S) is a base value dependent on the transformation used (S) and on the (base) sub-band considered 
foHnsertion, of the maximum weighting coefficients ensuring imperceptibility when a single coefficient of said 
sub-band is modulated, 
(3 is strictly greater than 2, and 

E[lwlP] is the mathematical expectation of the function IwlP 

56. Watermarking signal insertion method according to one of Claims 51 to 55, characterised in that said minimum 
value (a D ) is determined by using a calculation function also dependent on the variance (o x 2 ) of the coefficients 
to be modulated, the distribution of the watermarking signal (w) and possibly the variance (a n 2 ) of an additive noise 
decorrelated from the coefficients modelling any distortions applied to said coefficients. 

57. Watermarking signal insertion method according to Claim 56, characterised in that said calculation function is of 
the type 



a(Jo 2 x + o 2 n ) 
JbP- c 



where g x 2 is equal to the variance of the coefficients to be modulated, o n 2 is equal to the variance of the 
additive noiseand a, b and c are constants which depend on the distribution of the watermarking signal (w) and 
the required detection probability level (p). 

58. Watermarking signal insertion method according to one of Claims 51 to 57, characterised in that the watermarking 
signal (w) is a pseudo-random sequence with a predetermined distribution with a zero mean. 

59. Watermarking signal insertion method according to one of Claims 51 to 58, characterised in that it includes the 
following steps: 

- calculating (E602) a maximum value (a^N)) of the weighting coefficients for a length (P) of the watermarking 
signal equal to the cardinal number (N) of the set of modulatable coefficients (X) ensuring imperceptibility of 
the watermarking signal (w) at the limit of perceptibility; 

- calculating (E603) the length (P°) of the watermarking signal (w) for a minimal value (c^) of the weighting 
coefficient equal to said maximum value calculated (a v (N)) and for a predetermined detection probability level 
(P); and 

- comparing (E606) said calculated length (P°) with the cardinal number (N) of the set of modulatable coefficients 
(X). 

60. Watermarking signal insertion method according to Claim 59, characterised in that, if at the comparison step (E606) 
said calculated length (P°) is less than or equal to the cardinal number (N) of the set of modulatable coefficients 
(X), there is inserted a watermarking signal (w) of length (P) at least equal to said calculated length (P°) by mod- 
ulating a subset of coefficients with a cardinal number equal to said length (P) according to a linear model using 
a weighting coefficient (a) equal to said calculated maximum value (a v (N)). 

61 . Watermarking signal insertion method according to Claim 59, characterised in that, if at the comparison step (E606) 
the said length (P°) is greater than the cardinal number (N) of the set of modulatable coefficients (X), there is 
calculated the detection probability level (p') obtained for a minimum value (Oq) of the weighting coefficient equal 
to said calculated maximum value (a v (N)) for a length (P) of the watermarking signal (w) equal to the cardinal 
number (N) of the set of modulatable coefficients (X). 

62. Watermarking signal insertion method according to one of Claims 51 to 58, characterised in that it includes a step 
of calculating a threshold value (P*) of the length of the watermarking signal (w) determined so that said minimum 
value (oco) of the weighting coefficients is equal to said maximum value (a v ) of the weighting coefficient for a length 
(P) of the watermarking signal (w) equal to said threshold value (P*) and a predetermined detection probability 
level (p). 

63. Device for inserting a watermarking signal (w) in a set of coefficients representing a digital image (I), in which at 
least one subset of coefficients is modulated according to a linear model using a weighting coefficient (a), the 
weighting coefficient being greater than a minimum value (oq) determined as a function notably of the length (P) 
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pno e c no!T ark ' n9 Signa ' (W) 3nd 3 de,ec,ion Probability level (p), characterised in that it comprises means (500 
502, 503) of calculating a maximum value (a v ) of the weighting coefficient as a function of the length (P) of the 
watermarking signal ensuring imperceptibility of the watermarking signal (w). 

^o ,e ^ rk ' n9 S ' 9na ' inser1ion device according to Claim 63, characterised in that it also comprises means (500 
502, 503) of determining compatible values of the length (P) of the watermarking signal, of the weighting coefficient 
(a) and of the detection probability level (p) so that the value of the weighting coefficient determined (o) is less 
than or equal to sa.d maximum value (a v ) and greater than or equal to said minimum value (a D ) for said determined 
values of the length (P) of the watermarking signal and of the detection probability level (p). 

Watermarking signal insertion device according to one of Claims 63 or 64, characterised in that the means (500 
502, 503) of calculating said maximum value (a v ) use a function independent of said set of coefficients to be 
watermarked. 



^nnTno Lnaf ? "I™ 00 deV,Ce accordin 9 ,oone of Claims 63 *> 65, characterised in that it comprises means 
(500, 502, 503) for the spatio-f requency transformation of an image so that the coefficients are spatio-frequencv 
coefficients obtained by spatio-frequency transformation of an image, the calculation of said maximum value (a„) 
of the weighting coefficient depending on the type of transformation used (S). 

6? ' ^o e ITo, rkin9 Si9nal inSerti ° n dSViCe accordin 9 to Claim 66, characterised in that said calculation means (500 
502, 503) use a function of the type 

a v ( P, S, w) = — 

P 1/|J (E[\wf]) V * 



where 



ow^S) is a base value dependent on the transformation used S and on the (base) sub-band considered for 
the insertion, of the maximum weighting coefficient ensuring perceptibility when a single coefficient of the said 
sub-band is modulated, 
P is strictly greater than 2, and 

E[lw)P] is the mathematical expectation of the function livlP. 

68. Watermarking signal insertion device according to one of Claims 63 to 67, characterised in that said minimum 
value (ceo) is determined by using a calculation function also depending on the variance (of) of the coefficients 
to be modulated, the distribution of the watermarking signal (w) and possibly the variance (of) of an additive noise 
decorrelated from the coefficients modelling any distortions applied to said coefficients. 

69. Watermarking signal insertion device according to Claim 68, characterised in that said calculation function is of 
the type ' " ' 



a D( P) = W3ES 
J bP- c 



where of is equal to the variance of the coefficients to be modulated, of is equal to the variance of the 
additive noise and a, b and c are constants which depend on the distribution of the watermarking signal and the 
required detection probability level. 

70. Watermarking signal insertion device according to one of Claims 63 to 69, characterised in that the watermarking 
signal (w) is a pseudo-random sequence with a predetermined distribution with a zero mean. 

71 . Watermarking signal insertion device according to one of Claims 63 to 70, characterised in that it also comprises: 

- means (500. 502, 503) of calculating the length (P°) of the watermarking signal for a minimum value (ok) of 
the weighting coefficient equal to a maximum value (a v (N)) calculated for a length (P) of the watermarking 
signal equal to the cardinal number (N) of the set of modulatable coefficients ensuring imperceptibility of the 
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watermarking signal at the limit of perceptibility and for a predetermined detection probability level (p); and 
- means (500, 502, 503) of comparing said calculated length (P») with the cardinal number (N) of the set of 
modulatable coefficients (X). 

72 Watermarking signal insertion device according to Claim 71 , characterised in that it also comprises means (500. 
502 503) of inserting a watermarking signal (w) with a length (P) at least equal to said calculated length (P ). 
adapted to modulate a subset of coefficients with a cardinal number equal to said length (P) according to a linear 
model using a weighting coefficient (a) equal to said calculated maximum value (o v (N). if said calculated length 
(P°) is less than or equal to the cardinal number (N) of the set of modulatable coefficients (X). 

73 Watermarking signal insertion device according to Claim 71 , characterised in that it also comprises means (500 
" 502 503) of calculating a detection probability level (p') obtained for a minimum value (a D ) of the weighting coef- 
ficient equal to said maximum value (a v (N)) calculated for a length (P) of the watermarking signal equal to the 
cardinal number (N) of the set of modulatable coefficients. 

74 Watermarking signal insertion device according toone of Claims 63 to 70, characterised in that it comprises means 
(500 502 503) of calculating a threshold value (P*) of the length of the watermarking signal determined so that 
said min imum value (an) of the weighting coefficient is equal to said maximum value (a v ) of the weighting coefficient 
for a length (P) of the watermarking signal (w) equal to said threshold value (P*) and a predetermined detection 
probability level (p). 

75 Watermarking signal insertion device according to one of Claims 63 to 74, characterised in that it is incorporated 
. in a microprocessor (500), a read only memory (502) storing a program for inserting a watermarking signal in an 

; " image and a random access memory (503) containing registers adapted to store variables modified during the 
running of said program. 

76. Computer, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 . 

77. Computer, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22. 

78. Computer, characterised in that it comprises means adapted to implement the decoding method according to one 
of Claims 23 to 26. 

79. Computer, characterised in that it comprises means adapted to implement a method of processing an image ac- 
cording to one of Claims 47 to 48. 

80. Computer, characterised in that it comprises means adapted to implement a method of inserting a watermarking 
signal according to one of Claims 51 to 62. 

81. Digital image processing apparatus, characterised in that it comprises means adapted to implement a partitioning 
determination method according to one of Claims 1 to 21 . 

82. Digital image processing apparatus, characterised in that it comprises means adapted to implement the method 
of inserting a watermark according to Claims 22. 

83. Digital image processing apparatus, characterised in that it comprises means adapted to implement the decoding 
method according to one of Claims 23 to 26. 

84. Digital image processing apparatus, characterised in that il comprises means adapted to implement a method of 
processing an image according to one of Claims 47 or 48. 

85. Digital image processing apparatus, characterised in that it comprises means adapted to implement a method of 
inserting a watermarking signal according to one of Claims 51 to 62. 

86. Digital printer, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 . 
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95. Digital photographic apparatus, characterised in that it comprises means adapted to implement a method of , n 
sailing a watermarking signal according to one of Claims 51 to 62. implement a method of ,n- 

96 ' e^rt^^OTt^S*^1 > to " COmpriSeS means ada P ted t0 implement a partitioning determination method 

98 ' SS3 0^ 2 C 3t a o ra 26 enSed ^ th3t " COmPnSeS meanS ad3P,ed ,0 imp ' ement the dec ° d ^ method according to 
"• ^^T^^Tl " C ° mpriSeS m ~» 3dapted - -P'— a method of processing an image 
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'^a^n^or:^ 01 ' Char3CteriSed in that » is ^ implement a method or inserting a watermark 

103 'ofCrm: r 2 P 3 r To 26. ^ ChMCto ^ h that « is ^P«- «° implement a decoding method according to one 

""zssssz p c\2:^:r:r^ in that « is adap,ed to ,mp,em - a — °< p — 9 - image 

^.information storage means readable by a microprocessor, characterised in that i, stores a pro gra m imp.ementing 
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a partitioning determination method according to one of Claims 1 to 21 . 

107.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
a method of inserting a watermark according to Claim 22. 

1 Conformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
a decoding method according to one of Claims 23 to 26. 

1 (^Information storage means readable by a microprocessor characterised in that it stores a program implementing 
a method of processing an image according to one of Claims 47 or 48. 

110. Information storage means readable by a microprocessor characterised in that it stores a program implementing 
a method of inserting a watermarking signal according to one of Claims 51 to 62. 
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